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ABSTRACT 

We present a new submm/mm galaxy counterpart identification technique which 
builds on the use of Spitzer IRAC colors as discriminators between likely counter- 
parts and the general IRAC galaxy population. Using 102 radio- and SMA-confirmed 
counterparts to AzTEC sources across three fields (GOODS-N, GOODS-S, and COS- 
MOS), we develop a non-parametric IRAC color-color characteristic density distri- 
bution (CDD), which, when combined with positional uncertainty information via 
likelihood ratios, allows us to rank all potential IRAC counterparts around SMGs 
and calculate the significance of each ranking via the reliability factor. We report 
all robust and tentative radio counterparts to SMGs, the first such list available for 
AzTEC/COSMOS, as well as the highest ranked IRAC counterparts for all AzTEC 
SMGs in these fields as determined by our technique. We demonstrate that the tech- 
nique is free of radio bias and thus applicable regardless of radio detections. For ob- 
servations made with a moderate beamsize (~18"), this technique identifies ^85 per 
cent of SMG counterparts. For much larger beamsizes (>30"), we report identification 
rates of 33-49 per cent. Using simulations, we demonstrate that this technique is an 
improvement over using positional information alone for observations with facilities 
such as AzTEC on the LMT and SCUBA-2 on JCMT. 
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1 INTRODUCTION 

The discovery of a large population of bright sources at high 
redshift through the opening of the submillimetre (submm) 
and millimeter (mm) wavelength windows continues to have 
a profound impact on our understanding of galaxy evo- 
lution in the early Universe. These submm/mm-selected 
galaxies (hereafter SMGs) are characterized by large far- 
infrared luminosities (>10 12 £©), tremendous star forma- 
tion rates (> 500 AIq yr" 1 ), and a number density that 
is high compared to local ultra-luminous infrared galaxies 
(ULIRGs), indicating strong evolution within the popula- 



tion (e.g IScott et alj 1200 j ). SMGs are thought to repre- 
sent young, massive systems that exist in an epoch of rapid 
mass build-up and may be the predeces sors of modern-da y 
giant elliptical galaxies (see review by iBlain et al.ll2002l ). 
As a population, SMGs contribute a significant portion of 
the cosmic energy density at z—2-3 through some combi- 
nation of dusty sta r bursts and active galactic nuclei (AGN 



Hughes et al. 19981: [Alexander et al. 20051: Chapman et al 



20051: iKovacs et all l200rj: IPope et all l200rj: ICoppin et al 



20081 : Murphy et all 1201 ll ). Observations at submm/mm 
wavelengths benefit from a strong negative k-correction and 
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thus present a unique opportunity for an unbiased view 
of star formation out to extreme redshifts. Understanding 
SMGs is vital to any complete view of galaxy evolution. 



both ground-based surveys at 450-2000 /im (e. 


e. Smail et al. 


19971; Barger et all 19981; Hughes 


et all 1998 


; Cowie et al. 


20021; IScott et al.ll2002l; iBorvs et al 


ll2003l;lGreve et all 


20041; 


Laurent et al.l 120051; ICoppin et alT 


I2OO6I; fecott et al. 


2008; 


Perera et al. 2008: Austermann et al.ll2010l) and in balloon- 



and space-based surveys at 250 - 500 ixra llPascale et al.ll2008l : 
iDevlin et all 120091 : lEales et alj|2010l : lOliver et alj|2010l ). A 
complete understanding of the physical processes within 
these systems, however, requires that these observations be 
matched to multi-wavelength data, from radio to X-ray. 
Such multi-wavelength counterpart identification of SMGs 
is difficult due to the poor angular resolution (>10") of 
single-dish submm telescopes and the intrinsic faintness of 
their optical counterparts due to dust obscuration. A great 
deal of work has been invested into counterpart identifi- 
cation techniques, primarily using radi o continuum data 
from interferometric observations (e.g. Ivison et al] 
Chapman et all 120051; IPope et all | 2006l ; llvison et all 
Chapin et al]|2009l ;l Yun et al. I l2012l ~ Radio observatio 



2002; 



2007; 



particularly useful for counterpart identification due to 1) 
the low number density of radio sources, which minimizes 
the probability of a chance association, and 2) the tight 
correlation between radio and submm/mm emission (e.g. 
ICondon|[l992l : lYun. Reddv fc Condonll200ll ). both star for- 
mation tracers. Direct submm/mm interfero metric imaging 



of SMGs with high angular resolution (e.g. Younger et al.l 
l2007l . l2009l : ISmolcic et alll2012l : lBarger et al'llioij ) supports 
this radio-submm association. 

Unfortunately, the use of radio detections for counter- 
part identification has drawbacks. Due to the radio-submm 
correlation, radio counterparts will preferentially be found 
for the brightest SMGs, which may not be representative 
of the full SMG population. Interferometric observations 
have shown that a small percentage (~7-15%) of single ra- 
dio sources near SMGs are not actually associated with the 
submm emission. In cases where multiple radio detections 
are found within the submm beam, interferometry reveals 
that in ~80% of these instances, the submm emission is only 
associated with one of the radio sources (see Section 13.11 
Younger et ai]|2007, 2009; S molcic et al.ll2012l ; iBarger et all 



2012l L Finally, radio observations suffer from a strong pos- 



itive k-correction, which results in a bias against high red- 
shift galaxies, while submm emission remains equally de- 
tectable out to z~8 ijBlain et al.l [2002). Current radio ob- 
servations that overlap with submm surveys have placed the 



at ~40-100 per cent (llvison et all 20021. 120071: IWeifi et al. 


2009 


: IChaoin et al 


2009 


: iLindner et all 1201 ll: Yun et al. 


2012 


;ISmolcic et al. 


2012; Barger et al.l 2012f), depending on 



the depth of the radio observations. For example, given the 
shallowness of the available radio data in the COSMOS field 
(see Section [2}, we can predict that we will detect ~40 
per cent of SMGs in the radio based on submm counter- 
part studies of GOODS-N SMGs, which utilized deep radio 
data to build a cumulative distr ibution of radio counter- 
parts as a function of radio depth (Pope 2007). And, in fact, 
we find that 40 per cent of SMGs detected with single-dish 
observations in COSMOS are associated with radio detec- 



tions l|Aretxaga et al]|201ll . this work). This is further con- 
firmed by interferometric observations of AzTEC SMGs in 
COSMOS. lYounger et al] (|2007l . EH) imaged a sample of 
15 SMGs with the Submillimetre Array (SMA) and found 
that only 6/15 (40%) have a robust radio counterpart. Of 
the nine with no radio counterpart, seven are also unde- 
tected with MIPS 24 /im. On the other hand, 13/15 of these 
SMGs are detected with Spitzer Infrared Array Camera 
(IRAC; iFazio et al1l2004ar ). Of the two not detected, one 
is confused with a foreground source. Additionally, several 
studies have shown that SMGs are routinely detected in 



dlono et al.ll2006l; IWang et al. 20071. 2011 


; Hatsukade et al] 


2010; 


Tarmrra et al.l 


2010; Ikarashi et al. 


1201 ll). The high 



source dens i ty (^6 arcmin at the 1.4 /xJy level at 3.6 /im; 
IFazio et al] l2004bl ) of IRAC galaxies makes it impossible 
to use IRAC det ections directly for unique identifications 
ijPope et al.ll2006h . however, the selection of IRAC counter- 
parts is expected to suffer less bias against high redshift 
sources and may be more representative of faint SMGs (see 
Section [572]) . 

More recently, two studies have examined follow-up in- 
terferometry of SM Gs detected with LABOCA (870 /ira; 
' Siringo et all [2009) and SCUBA/SCUBA-2 (850 /im)- 
Smolcic et al.l |2012l) analysed a set of 34 LABOCA sources 



in COSMOS with interferometry from various sources, in- 
cluding 26 new detections from the IRAM Plateau de 
Bure Interferometer (PdBI). Using a radio catalog extended 
to include lower signal-to- noise (S/N) sources, they found 
that 18/34 (~50%) SMGs were detected at 20cm and the 
same fraction, though a different subset, were detected with 
IRAC. This study also determines redshifts to their sources 
and concludes that the redshift distribution of SMGs may 
peak higher than previously thought, including a large tail 
of sources at very high redshift (z>3) that will be difficult 
to detect in even the deepest radio surveys due to the k- 
correction. It should be noted, however, that this study did 
not require the LABOCA SMGs to be detected at high sig- 
nificance. If we limit the sample to high significance inter- 
ferometric detections (>4.5<r), then 10/16 (62%) SMGs are 
detected in IRAC, while, if we look only at LABOCA sources 
which are also detected wit h AzTEC, 8/1 (80% ) have IRAC 
counterparts. Additionally. IBarger et all (|2012T l looked at a 
small, but homogenous sample of 12 SMGs in GOODS-N 
and found 16 independent, high significance detections with 
the SMA. Of these 16, all are detected in an ultra-deep radio 
catalog and 15 are detected in the IRAC bands. 

In this study, we develop a new technique for counter- 
part identification which takes advantage of the high sensi- 
tivity of Spitzer IRAC and does not rely on radio or MIPS 
24/im observations. Our starting point is the IRAC color- 
color space diagram first developed by lLacv et al] (|2004l ) 
for identifying AGN. IRAC colors, and this diagram in par- 
ticular, have been shown to also be useful for SMG coun- 
terpart identification, with <20 per cent of SMGs overlap- 
ping with AGN in IRAC color- color space (|Pope et al]|2006l ; 
iBiggs et al] l201ll ; lYun et all 120081 . I2012T ). We take these 
studies a step farther and, using a large sample of SMGs with 
secure radio or SMA counterparts, develop an IRAC color- 
color density distribution. This multi-dimensional color 
prior is combined with posi tional uncertainty informatio n 
via the likelihood ratio (LR; ISutherland fc Saunderslll992l ). 
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providing a method for ranking potential IRAC counterparts 
around SMGs. Previous studies utilizing likelihood ratios 
for counterpart identification can be found in the literature 
across all wavelength regimes; for example, IChapin et al.l 
Hon]) combined priors based on radio, MIPS, and IRAC 
fluxes and colors with positional information to identify 
count erpar ts to BLAST and LABOCA sources. ISmith et al.l 
l|201lf l and iKim et all (|2012t ) similarly used individual flux 
and color priors to identify counterparts to Herschel SPIRE 
sources. The methods used in these studies necessarily as- 
sume that the various priors u sed are independent , though 
some or all may be correlated (|Chapin et al.ll2~01ll ). In this 
work, we incorporate correlations between different flux 
bands and colors by using a characteristic density distri- 
bution (CDD) rather than individual priors. 

Since our method only requires Spitzer IRAC obser- 
vations, it can be applied in numerous fields that lack deep 
radio observations and is not as strongly biased against high 
redshift sources as techniques which rely on radio data. 
Many of the fields in which deep IRAC data already ex- 
ists will be prime targets for future submm surveys and we 
show that the technique developed here is useful not only 
for existing submm surveys, but also for upcoming surveys 
with the LMT and SCUBA-2. Though the particular char- 
acteristic density distribution developed in this study is only 
applicable to source populations with a redshift distribution 
and spectral energy distribution (SED) similar to sources 
selected at 850-1200 /im (see Section l3.2.ip . we posit that 
a modified version of this technique can be developed for 
Herschel SPIRE sources, which suffer from similar chal- 
lenges due to low resolution. 

This paper is organized as follows. In Section [2] we de- 
scribe the data used in this study. In Section [31 we develop 
the counterpart identification technique, which we use in 
Section [4] to calculate an LR and its corresponding relia- 
bility for all IRAC counterpart candidates across the three 
fields, including an analysis of the resulting ranking via a 
blind-test using our radio- and SMA-identified counterpart 
sample. In Section [S] we discuss our technique in the context 
of future submm surveys and address the issue of possible 
bias introduced through the use of radio counterparts in 
developing the IRAC color-color characteristic density dis- 
tribution. Section [6] contains our conclusions. 



2 DATA 

To develop our new technique for selecting IRAC counter- 
parts to SMGs, we utilize a statistically signific ant sample of 
272 S MGs detected at 1.1 mm with AzTEC (|Wilson et all 
2008) from three fields: the Cosmic Evolution Survey (COS- 
MOS), the Great Observatories Origins Deep Survey-North 
(GOODS-N), and -South (GOODS-S). These three fields 
have extensive multi-wavelength data ranging from X-ray 
to radio, much of which was observed at high resolution. We 
take advantage of observations from the Very Large Array 
(VLA) at 1.4 GHz and interferometric observations taken 
with the SMA to identify secure counterparts for a subset 
of the AzTEC SMGs; the IRAC counterparts are then used 
to derive a characteristic density distribution utilizing the 
four Spitzer IRAC bands (3.6, 4.5, 5.8, and 8.0 pm), which 



is described in the next sections. We describe each data set 
below. 

AzTEC: The AzTEC observat ions used h er e are 
described in th e following papers : IPerera et all (|2008h 
for GOODS-N IScott , et""aU (120101) for G OODS-S, and 
IScott et alJ (|2008l ) and lAretxaga et ail (|201lf ) for COSMOS. 
Data reductio n was performed using the standard A zTEC 
pipeline (see IScott et all |200S|; IDownes et alJ I2012T I. The 
reader is referred to IDownes et al.l ( 20121 ) for the most up- 
dated source lists for the GOODS fields. Observations in 
GOODS-N were taken on the JCMT with uniform cover- 
age to a depth of ~1.3 mjy beam -1 and a full width at 
half max (FWHM) of 18". GOODS-S and COSMOS were 
imaged using the Atacama subm illimetre Telescope Experi- 
ment (ASTE: lEzawa et alj2004 ) to a depth of ~0.6 and ~1.3 
mjy beam -1 , respectively. These observations were done at 
a lower angular resolution and we measure a FWHM of 30" 
for GOODS-S and 33" for COSMOS by fitting the post- 
filtered PSFs of each set of observations with a Gaussian. 
Our sample contains 36, 47, and 189 SMGs in GOODS-N, 
GOODS-S, and COSMOS, for a total of 272 sources. 

VLA: Radio observations at 20 cm (1.4 GHz) were ob- 
tained for all fields using the VLA with an angular resolution 
of ^3.5". Our deepest radio catalog is in GOODS-N, with 
a sensiti vity of 3.9 fj,Jy beam - 1 (la) and 1,230 5a radio de- 
tections (|Morrison et alj|20ld ). Our GOODS-S catalog was 
made using the "Search and Destroy" (SAD) procedure in 
AIPS to identify all sources greater than 3a, which were 
then cross-matched with the IRAC SIMPLE and GO ODS 
catalogs to exclude noise peaks (see lYun et "ail 12012ft . The 
GOODS-S catalog has an rms sensitivity of 8.5 /xJy beam -1 
and ~ 1 ,500 sources. For COSM OS, we use the >4.5a cat- 
alog of ISchinnerer et al. 

I (|2010r i. which has ~3,000 sources 
to a la sensitivity of 10 /xjy beam - . We demonstrate in 
Section l5.3l that the sensitivity of the radio catalog has little 
impact on our results. 

SMA : Fifteen AzTEC/ JCMT sources in the COS- 
MOS field w ere the targets of two foll ow-up observations 
by the SMA jYounger et alj|2007l . l2009h . All were detected 
at 890 fim with positional accuracies of ~0.2". 6/15 sources 
are detected in the radio, while nine are not. We include 
the three radio-undetected sources that were detected in all 
four IRAC bands in our analysis and in our derivation of 
the IRAC color-color CDD below. 

Spitzer IRAC : GOODS-N, GOODS-S, and COSMOS 
were imaged with all four IRAC bands with an angular reso- 
lution of 2" and a positi onal accuracy of ^0. 2". In GOODS- 
S, the SIMPLE survey l|Damen et all 120111 ) covers ~1,600 
arcmin -2 to a 5a depth of 1.0, 1.2, 6.0, and 6.7 /iJy for 
bands 3.6, 4.5, 5.8, and 8.0 /im. The SIMPLE IRAC catalog 
comprises all 3.6 fim detections with S/N > 5. Deeper obser- 
vations were done wit h a 5a sensitivity of 0.1, 0.2, 1.5, and 
1.6 ^Jy in GOODS-N (|Dickinsodl2003l ; iTreister et alj|2006ft 
and shallower observations were take n with depths of 0.9 , 
1.7, 11.3, and 14.6 ^Jy in S-COSMOS (|Sanders et alj|2007ft . 
with a catalog containing all sources detected at 3.6 fim at 
the >1 /j,Jy level. Based on studies in the GOODS fields, 
these sensitivities allow for the detection of galaxies with 
stellar masses greater than 10 10 Mq a nd a typical age of a 
few hundred million years up to 2 ~ 6 (|Yan et alJl2"b"oa ). 
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3 DERIVING THE CHARACTERISTIC 
DENSITY DISTRIBUTION 

Our submm galaxy counterpart identification technique re- 
lies on two pieces of information: the positional uncertainty 
of a given AzTEC detection and the general distribution of 
IRAC colors for the SMG population. The latter is obtained 
from our 'training set', the largest sample of radio- and 
SMA-identified counterparts to SMGs selected at 1.1mm 
used to date in a counterpart identification study. We de- 
scribe the selection of secure radio counterparts and the de- 
velopment of the IRAC color-color CDD of SMGs in this 
section. 



3.1 Identifying Radio Counterparts to SMGs: the 
Training Set 

Due to the typically low S/N (<10) and low angular resolu- 
tion of submm/mm observations with single-dish telescopes, 
the positional offsets between the detected emission and the 
true position of the source are large due to noise fluctuations 
in the the submm maps. We can characterize these posi- 
tional uncertainties though simulations wherein simulated 
sources of known flux densities are inserted into the map 
one at a time at random positions. From this we build a 
cumulative positional probability distribution P(>D; S/N), 
which is the probability that a source will be detected at a 
distance greater than D from its true position as a func- 
tion of S/N, where the noise is the instrument noise of 
the observation (see iPerera et al.l 120081 ; IScott et all l20ld ; 
lAretxaga et al.l 1201 if ). The distributions determined from 
these simulations are co nsistent with the analytical solution 
from llvison et alj ((2007), where 

P(> D; S/N) = j D 1 exp (j^rj dD ' M 
with 

0.6FWHM 

^ S/N (2) 

Using this formula, we can identify all potential ra- 
dio counterparts around each SMG using a variable search 
radius based on the S/N of the AzTEC detection. To be 
conservative, we select all radio sources within a search 
radius where P(>D; S/N) goes to one per cent. For the 
lowest S/N sources considered here (S/N = 3.5), this cor- 
responds to a search radius of 9.5", 16", and 17.6" for 
SMGs in GOODS-N, GOODS-S, and COSMOS, respec- 
tively. Any AzTEC source that has a single radio source 
within the appropriate search radius is included in our 
training set. From interferometric follow-ups to AzTEC and 
LABOCA SMGs, we expect that the false association rate 
for single radio sources detected within the submm beam 
to b e ~7-15% ((Younger et all I20071 . 120091 : ISmolcic et ail 
l2012h . Using a variable search radius based on the posi- 
tional uncertainties, as in this work, should decrease this 
rate. An analysis of the twelve radio-detected SMGs in 
our training set that have interferometric data available 
(AzTEC/GN3, AzTEC/GN5, AzTEC/GN14, AzTEC/C3, 
AzTEC/C12, AzTEC/C13, AzTEC/C14, AzTEC/C18, 
AzTEC/C22, AzTEC/C38, AzTEC/C98, AzTEC/C145) 



shows that we have selected the correct radio counter- 
part in all cases except AzTEC/GN14 (see Table |A"T1 
Younger et ai1l2007l . 120091 ; ISmolcic et aT]l2012l : iBarger et all 
120121 ). AzTEC/GN14 has a faint radio counterpart that is 
below the detection limits of our catalog, in addition to an 
unrelated radio source within its search radius. 

Traditionally, the robustn ess of a radio count erpart is 
determined by the P-statistic ((Downes et al.l fl986T ) . with a 
P-stat < 0.05 considered a robust counterpart and 0.05 < 
P-stat < 0.2 considered a tentative counterpart. We identify 
a total of 107 isolated radio counterparts to our SMGs, 18 in 
GOODS-N, 15 in GOODS-S, and 74 in COSMOS. Of these, 
95 (79%) have P-stat < 0.05, while 22 (21%) have 0.05 < P- 
stat < 0.2. We elect to include these tentative counterparts 
in our training set and have found that they do not signifi- 
cantly influence the resulting IRAC color distribution. Our 
ra dio counterpa r ts are in good agreement with those found 
in lChapin etaD (|2009l ) for GOODS-N and lYun et ail (|2012f ) 
for GOODS-S. This work presents the formal analysis of the 
ra dio counterparts to A zTEC/COSMOS sources described 
in lAretxaga et al.l (120111 ) (Tables |ATIA"2)) . 

IRAC counterparts are obtained by searching within 
2" of the radio positions, which is possible due to the ex- 
cellent positional accu racy of the VLA observations (<l"; 
iKellerman et al.l 120081 ') . For determining the IRAC color- 
color CDD, we require that IRAC sources be detected in all 
four bands. This leaves 96/107 (90%) radio-detected SMGs 
with complete IRAC data. Of the remainder, ten are com- 
pletely undetected and one is detected in only 3 bands. 
To this sample, we add an additional six SMGs with se- 
cure IRAC counterparts (detected in all four bands) deter- 
mined from the SMA imaging that are radio-undetected, 
have multiple radio counterparts, or are not in the current 
AzTEC/COSMOS catalog. This gives a total of 102 IRAC 
counterparts which we use to build an IRAC color space 
CDD. The list of radio and IRAC counterparts to radio- 
detected AzTEC sources! can be seen in Table EU includ- 
ing the search radii used and associated P-statistics. Those 
AzTEC sources not included in the training set due to lack 
of an IRAC counterpart are indicated. 

We find that 8 (22%), 4 (9%), and 15 (8%) SMGs in 
GOODS-N, GOODS-S, and COSMOS have multiple radio 
counterparts 1 within their search radii (Table |A"2)) . SMGs 
with multiple radio counterparts are not included in the 
training set due to ambiguity in the source of the mm emis- 
sion. 

3.2 IRAC Color-Color Space 

Follo wing the lead of lLacv et all (|2004h and lYun et ail (|2008l . 
l2012h , we create an IRAC color-color diagram with the col- 
ors log (§f-f ) versus log (fM-J . The 102 IRAC counterparts 
in our training set can be seen in IRAC color-color space in 
Figure[T] along with a random sample of 10,000 IRAC galax- 
ies from our three IRAC cat alogs. The dotted lines represent 
the color cuts presented in lYun et al.l l|2008l ). which encase 
the bulk of IRAC colors for securely identified SMGs. 

1 A list of AzTEC sources and radio/IRAC coun- 
terparts, including photometry, is available at 
www. astro .umass . edu/aztec/CounterpartlD/counterpartid.html 
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Figure 1. The radio- and SMA- identified IRAC counterparts to 
SMGs that are used as a training set in this study. In addition to 
the 102 IRAC counterparts detected in all four IRAC channels, 
4 counterparts that are missing one or more IRAC detections are 
plotted as upper limits. These sources are not part of the train- 
ing set. The dotted lines denote the IRAC color cuts suggested 
for SMC counterpart identification in lYun et al.l l|2008l ). The gray 
dots represent a random subsamplc of ~10,000 IRAC field galax- 
ies. 



The following describes the building of the IRAC color- 
color CDD for SMGs, starting with a brief discussion of the 
physical mechanisms responsible for the IRAC color-color 
properties of SMGs. 



3.2.1 The Physical Mechanism Behind the IRAC 
Color-Color Distribution of SMGs 

Due to the large redshift range at which SMGs have been 
detected thus far (z~l-6 with a median redshift of ~2-2.5; 



e.g. Smail et al.l2004; Chai 


smart et al. 2005 


: Aretxaga et al. 


20071; IChapin et alJ 120091; 


Aretxaea et al. 


20 111; Yun et al. 



20121 ; ISmolcic et al.ll2012f ). the IRAC bands probe a handful 
of spectral features over the rest-frame wavelength range of 
0.5 to 4.0 /im. In a typical star forming galaxy, this range in- 
cludes the stellar photospheric feature at rest-frame 1.6 /im 
(the "stellar bump"), which is the continuum emission from 
the old, low mass stars that form the bulk of the stellar 
mass of the galaxy, and, at wavelengths >3 /im, emission 
from polycyclic aromatic hydrocarbons (PAHs). In the case 
of extremely young galaxies, emission from OB stars can 
shift the bulk of the stellar emission to shorter wavelengths, 
depending on dust content. Systems that are dominated by 
AGN show power-law emission from hot dust in this wave- 
length regime. 

To determine which of these features predominantly 
give rise to the observed IRAC colors of SMGs, we gener- 
ated SMG galaxy mode ls using the GRASIL SED modeling 
code (Silva et al. 1 1 19981 ) based on the following properties: 
SMGs are dominated by stellar emission as opposed to dust 
heated by AGN in the majority of cases ([A lexande r et "ah! 



is ~40-100 Myr l|Greve et afll2005l ; iNaravanan et al]|20ich 

and submm galaxies are dusty systems (Blainet al. 2002) 
which contain a pre-exi s ting stellar population |Smail et al 



2005; Pope et al. 



20081; iMenendez -Delmestre et al.l 120091 : 



Johnson et al the typical timescale of the SMG phase 



2004 1; iBorvs et al.ll2005l ; iDve et al l [20081 ; iMichalowski et al 
20121 ). We create several GRASIL models with galaxy ages 
ranging from 1-2 Gyr, each undergoing an exponentially de- 
clining starburst with timescales ranging from 10-100 Myr 
and e- folding times ranging from 20-100 Myr. Each model 
contains a moderate amount of dust with a dust to gas ra- 
tio of 5 per cent and an optical depth at 1 /jm of 0.01. For 
extremely youn g star forming syst ems, w e com piled a set 
of models from lBruzual fc Chariot! (BC03; 120031 ) with ages 
1-100 Myr, both with no ext inction and wi th Av = 4 using 
the Calzetti extinction law (|Calzettill200ll ). Figure [2] (top) 
shows GRASIL models with a galaxy age of 1 Gyr with star- 
burst ages of 10, 20, 50, and 100 Myr as well as two 1 Myr 
BC03 model SEDs, with and without dust extinction. The 
BC03 1 Myr model without extinction is only for reference 
as all SMGs are dusty systems. When extinction is added 
to extremely young systems, the emission from young stars 
is attenuated, creating a pseudo-stellar bump that is shifted 
to shorter wavelengths. In both model sets, the stellar bump 
is prominent by ~10 Myr, making it the dominant feature 
in the IRAC bands at the redshifts typical of SMGs. 

Of course, the stellar bump is not unique to SMGs; 
previous studies have established that the IRAC bands are 
a useful redshift selection tool for star forming galaxies at 
z>\ due to the perva siveness and robustness of the stellar 
bump and its shape llPope et alJ 12006) ; iFarrah et alJ [20081 : 
lLonsdale etafl 12009] ; ISorba fc Sawickif l201Ch . There is an 
additional complication in that at high redshift (z>4), the 
IRAC bands begin to sample only one side of the stellar 
bump, which can mimic a power-law shape, thus making 
it difficult to determine if these sources are star formation- 
dominated or AGN-dominated. We argue that AGN are rare 
and it is already known that there is some overlap in SMG 
and A GN color in IRAC color-color space (|Yun et al.ll2008l . 
120121). For further disc ussion on the IRAC colors of AGN, 
see lDonlev et al.1 i|2012T ). 

The redshifted color-color tracks of our GRASIL models 
are shown in Figure [2] (bottom) overlaid on the IRAC colors 
of our SMG training sample. For comparison, the IRAC col- 
ors of a random sample of field galaxies are shown as well. 
It can be seen from the plots in Figure [5]that the separation 
in IRAC color-color space between SMGs and the bulk of 
field galaxies is predominantly a redshift effect. As such, we 
are essentially using the IRAC colors to separate the high 
redshift SMG counterparts from the unrelated low redshift 
population. 



3.2.2 Kernel Density Estimation 

Our characteristic density distribution is developed using 
kernel density estimation, which provides a non-parametric 
estimator of t he underlying probability density function 
of a data set (iParzenl 1 19621 ; iPisanil 1 19931 ; IVio et al.l Il994l ; 
iMerritt et al.lll994l '). Functionally, this process convolves the 
IRAC colors of our training set with an optimized kernel 
to create a smoothed distribution which approximates the 
underlying parent distribution. Cross-validation is used to 
optimized the kernel, with the optimal kernel width being a 
function of both the sample size of the data and the broad- 
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Figure 2. (Top) Starburst models of various ages with the IRAC 
filters at 2=2 overlaid for reference. The GRASIL models repre- 
sent a typical SMG with an underlying 1 Gyr old population 
and a starburst with ages of 10-100 Myr; these models show the 
ubiquitous stellar bump and a PAH feature at 3.4 /im. The BC03 
models represent an extremely young lMyr starburst. An extreme 
starburst with no dust (A„ = 0) is not physical for SMGs, which 
are always dusty systems, so it is only shown for reference. When 
dust is added (A„ = 4), the emission from young stars is attenu- 
ated, forming a pseudo-stellar bump even in these extreme young 
starbursts. In all models we examined, the stellar bump is evi- 
dent by lOMyr. (Bottom) The redshifted IRAC color-color tracks 
of the GRASIL 50Myr SMG model between z=l-6 overlaid on 
the IRAC colors of the training set and a random subsample of 
IRAC field galaxies. 



ness (i.e. the standard deviation) of the intrinsic distribu- 
tion. To build the CDD, we use the 102 radio- and SMA- 
identified IRAC counterparts of our training set described 
in Section T3. 2 1 The steps for cross-validation are as follows: 

(i) Divide the training set in half randomly, forming a 
'working' set and a 'cross-validation' set in IRAC color-color 
space. 

(ii) Incorporate the photometric uncertainties into the 
working set by replacing each datapoint within IRAC color 
space with a normalized 2-d Gaussian with a standard de- 
viation defined by the lcr uncertainties on the colors. This 
can be seen in Figure [3] overlaid with the datapoints that 
form the cross-validation set. 



log(S 5 . 8 /S 3 . 6 ) 



Figure 3. An example of one cross-validation realization. The 
white contours represent the "working" set, a randomly chosen 
half of the training set. The datapoints in the working set have 
been replaced with normalized 2-d Gaussians in order to incor- 
porate their photometric uncertainties. The red pluses represent 
the "cross-validation" set. 



(iii) Next, select a kernel, which of its parameters to vary, 
and a suitable parameter space. We choose a 2-d Gaus- 
sian kernel and allow its standard deviations, a i s 5 s \ and 

og \ . B ) 

a i s s o \ , to vary independently. Preliminary tests showed 

° S \ s 4.5 ) 

that the optimized kernel has negligible rotation and so we 
do not allow the kernel to rotate in the following steps for 
computational simplicity. Convolve the working set with the 
kernel while varying its parameters over the parameter space 
chosen. 

(iv) For each set of parameters, compare the convolved 
working set to the cross-validation set. The optimized kernel 
parameters are those which maximize the likelihood that 
both the convolved working set and the cross-validation set 
are drawn from the same parent distribution. This likelihood 
is calculated as 



L = liiC I log 



where 



S5.8 
S3.6 



,log 



Ss.o 

S4.5 



is the value of the convolved working set at the ith datapoint 
in the cross-validation set. 

(v) Finally, repeat steps i-iv to build a representative dis- 
tribution of optimized kernel parameters. We perform the 
cross-validation procedure 500 times and find that the distri- 
bution of optimized Gaussian widths does not significantly 
vary in the range of 0.07 to 0.21 for both colors, which indi- 
cates that the final kernel estimation results are not overly 
sensitive to the exact size of the Gaussian kernel. 



We adopt the mean values a 



^ 1 . r» / 



(3) 



, /s 58 \=0.14 and 



f, /s sn \=0.15 for the final convolution of the 102 IRAC 



counterparts in the training set with an optimized kernel. 
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Our final IRAC color-color CDDsmcE 

for submillimetre 

galaxies can be seen in greyscale in Figure [4] overlaid with 
the original IRAC datapoints used to create it. 

To use this characteristic density distribution to identify 
counterparts to the rest of our SMGs, we must also deter- 
mine the IRAC color distribution for the full IRAC popu- 
lation. We use the same kernel density estimation discussed 
above on all IRAC galaxies that are detected in all four 
IRAC bands across the three fields, ~295,000 sources, and 
find a smaller range in optimal kernel widths, 0.02-0.06. We 
adopt the mean value of 0.04 for both standard deviations 
of the Gaussian kernel. The final distribution CDDirac, 
shown in red contours in Figure 21 peaks in a bluer regime 
of IRAC color-color space than the distribution for SMGs. 
This is expected since the full IRAC sample is dominated 
by z<l galaxies. At these lower redshifts, the IRAC bands 
sample the more complex PAH emission features, and the 
variety in observed galaxy SEDs at low redshifts accounts 
for the large scatter in IRAC color-color space of the general 
IRAC population. Regardless, the distribution of the field 
IRAC galaxies is distinct from that of the SMG population, 
which allows us to make use of this method for counterpart 
identification. 



4 IDENTIFYING IRAC COUNTERPARTS TO 
SMGS 

4.1 The Likelihood Ratio and Reliability Factor 

The relative probability that an IRAC galaxy is the coun- 
terpart of an SMG as opposed to being an unrelated fore- 

C round/background source is given by the likelihood ratio 
Sutherland fc Saunders! fl992l ): 



A = 



CDDsmg (log log (ff)) 

CDDxhao (l°g (ft), log fe)) 



P(D; S/N) (4) 



where CDD SM G (log (|M) , log (fff)) is the value 
of the CDD of SMGs at the location of a po- 
tential IRAC counterpart in IRAC color-color space, 

CDDirac (log (§|<f ) , log (f*f )) is the value of the CDD 
of the general IRAC population for that potential counter- 
part, and P(D; S/N) is given by Eqn. [T] integrated over D-e 
to D+e, where D is the distance of a given IRAC galaxy from 
the centre of the detected submm emission and e=0.2" is the 
positional accuracy of the IRAC observations. P(D; S/N) 
represents the probability that the true SMG position is 
found at distance D, based on the S/N and FWHM of the 
submm detection. 

The likelihood ratio provides a means of ranking IRAC 
galaxies in the proximity of SMGs based on their colors and 
positions. Since there is no IRAC color combinations that is 
completely unique to SMGs, this ranking is relative within 
a set of potential counterparts. To determine the statistical 
significance of any given likelihood ratio, we calculate its la 
uncertainty by including the measured errors in the IRAC 
colors. The ratio CDDsmg /CDDirac can be seen in Fig- 
ure [S] in black, dash-dot contours, overlaid on CDDsmg in 



greyscale and CDDirac in red, solid contours. Due to the 
rapid falling off of both the SMG and IRAC CDDs in IRAC 
color-color space, a floor value of log / = -5.7 was imposed on 
both distributions when plotting the ratio to suppress noise 
features in areas with sparse data. This floor was not used 
in the calculation of the likelihood ratio and we have verified 
that these noise features are not affecting our counterpart 
identification results. The peak of the CDDsmg /CDDirac 
contours show the area in which the color prior portion of 
the likelihood ratio is at a maximum. 

Addit ionally, we adopt the r eliabi lity factor, as de- 
scribed in ISutherland fc Saunders] (|l992T l. which provides a 
single statistic for a given potential counterpart which in- 
corporates its individual LR, the LRs of all other potential 
counterparts around that SMG, and the probability that 
the true counterpart is detected. The reliability factor for 
potential IRAC counterpart j takes the form 



A, 



E,A l + (l-A) 



(5) 



where A; is the sum over the likelihood ratios of all po- 
tential IRAC counterparts around a given SMG and X is 
the probability that the true IRAC counterpart is included 
in the four-band IRAC catalog we are using. Though in- 
terferometric studies indicate a high probability that the 
majority of SMGs are detected a t ^1 ujy at 3.6 and 
4.5 fim (e.g. 



l2010l : llkarashi et al.ll201lD . our technique requires detections 
at 5.8 and 8.0 /im as well. These bands are less sensitive than 
the two shorter IRAC channels, though our training set sug- 
gests a high percentage (99%) of counterparts are detected 
in all bands. The overall fraction of IRAC sources not de- 
tected in all bands therefore provides an estimate of X. 

To control for possible biases in the properties of IRAC 
galaxies that are nearby SMGs as opposed to the general 
field population, we elect to estimate X only from IRAC 
galaxies within the vicinity of SMGs, including both our 
training sample and the radio-undetected SMGs. To do this, 
we identify all IRAC sources within the search radii de- 
scribed in Section 13.11 and listed in Tables IAUA3I We find 
that the fraction of IRAC sources detected in all four bands 
is 0.83, 0.78, and 0.73 for GOODS-N, GOODS-S, and COS- 
MOS, respectively, and we use these as approximations for 
X in Eqn [5] We note that this fraction only changes by <10 
per cent if we include all IRAC catalog sources, suggest- 
ing no strong biases to this fraction from redshift selection 
effects. 

In principle, Xi for a given SMG in Eqn.[5]runs over 
the entire IRAC four-band catalog. In practice, however, 
Ai decreases rapidly with increasing distance between the 
mm peak emission and the IRAC position (Eqn [4}. For this 
reason, we limit the sum in the denominator in Eqn to 
IRAC galaxies within the search radii for the SMG, defined 
where P(>D; S/N) — 0.01, to simplify computations. 



ivHjS are detected a t -gi ujy at 3.o and 
Younger et all 120071 . 120091 ; lliatsukade et al.1 



4.2 Counterpart Identification via the Reliability 
Factor 

4-2.1 Verification of the Method: the Training Set 

We evaluate our technique with a blind-test against our 



2 Availableatww.astro.umass.edu/aztec/CounterpartID/counterpaiBak4ioitfllid SMA-identified training set. All IRAC galaxies 
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Figure 4. The final IRAC color-color characteristic density distribution for SMGs (CDDsmg\ white image contours) overlaid with the 
radio- and SMA-identified IRAC counterparts used to create it, as in Figure [l] Additionally, the CDD of all IRAC galaxies is shown 
(CDDirac'i re d contours). 



Table 1. Statistical results of counterpart identification. Columns 2-5 show the 
distribution of ranks for the training set per field. Being ranked first is considered 
a correct identification. Additionally shown are the number of robust and tentative 
reliabilities in the training set and the radio-undetected set. 



Field 


Ranked 
1st 


Ranked 
2nd 


Ranked 
3rd 


Ranked 
3rd+ 


R^0.8 (R^0.6) 
Training Set a Radio-Undetected 


GOODS-N 
GOODS-S 
COSMOS 


11 (85%) 
5 (33%) 
33 (49%) 


2 (15%) 
5 (33%) 
15 (22%) 




2 (13%) 
8 (11%) 




3 (20%) 
12 (17%) 


6/6 (9/9) 2 (8) 
2/3 (3/7) 6 (14) 
11/12 (18/28) 21 (48) 



a Thc number of correctly identified counterparts in the training set with a reliability factor 
above 0.8 (0.6) compared to the total number of IRAC galaxies assigned a reliability above 
that value. 



within the appropriate search radius of each SMG are ranked 
using the likelihood ratio (Eqn. 2| and assigned a reliability 
(Eqn. [5]). The likelihood ratio, its la errors, and reliabil- 
ity of the radio- or SMA-identified counterpart as well as 
additional top ranked potential IRAC counterparts for the 
training set can be seen in Table \KT\ 

In GOODS-N, 11/13 (85%) IRAC counterparts are 



identified correctly, i.e. ranked first among the other pos- 
sibilities using the likelihood ratio. Three of these sources 
(AzTEC/GN3, AzTEC/ G5, and AzTEC/G N14) are con- 
firmed by interferometry (jBarger et alj|2012h . In GOODS- 
S, 5/15 (33%), and in COSMOS, 33/68 (49%), are 
identified correctly, including one radio-undetected SMG 
(AzTEC/C4), one multiple radio system (AzTEC/C42, see 
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Figure 5. The ratio C DD smg /CDDjrac in black, dash-dot 
contours, overlaid on CDDsmg m the white image contours, and 
CDDjuac in red, solid contours. The datapoint shown here are 
the IRAC counterparts in the training set which were identified 
correctly by their likelihood ratio. 



below), and two radio-detected SMGs (AzTEC/C18 and 
AzTEC/C22) which were imaged with the SMA. The cor- 
rectly identi fied counterparts, as well as SMA-confirmed 
counterparts (| Younger et al.ll2007l . l2009l ). can be seen in Fig- 
ure [SI overlaid on contours showing CDDsmg, CDDirac, 
and the ratio CDDsmg /CDDirac ■ The difference in the 
success rates of GOODS-N versus COSMOS and GOODS- 
S stems from the difference in beamsizes: 18" in GOODS- 
N compared to 33" and 30" in COSMOS and GOODS-S. 
Possible explanations for the more subtle differences be- 
tween COSMOS and GOODS-S are discussed in Section 15751 
The distribution of counterpart rankings (1st, 2nd, 3rd, and 
3rd+) in each field can be seen in Table [1] 

The distribution of reliabilities for the radio- and SMA- 
identified training set can be seen in Figure [S] A large R 
value indicates a high probability that a given IRAC galaxy 
is the true counterpart as well as a high probability that only 
a single galaxy is responsible for the submm/mm emission. 
A low R value indicates a low probability of being a counter- 
part and/or that the corresponding SMG consists of multiple 
galaxies blended within the beam (see Section l5.1|) . Figure[B] 
shows that for R^0.5, all radio- or SMA-identified IRAC 
counterparts are identified correctly, whereas for R^0.5, 
many radio-identified IRAC counterparts are not. We sug- 
gest a conservative cut-off in the reliability factor would be 
at R^0.8, where there are only two misidentifications out of 
nineteen. This cut-off is the same as th at found in simila r 
studies utilizing the reliability factor (e.g. lSmith et al.ll2011f) . 
Additionally, an analysis of the field-by-field reliability fac- 
tors suggests that extending the cutoff to R^0.6 is appro- 
priate in fields analogous to GOODS-N, where 9/9 (100%) 
IRAC galaxies assigned R^0.6 are the correct counterpart. 
The range 0.6^ R ^0.8 is more appropriately named a 'ten- 
tative' counterpart list in COSMOS with only 18/28 (64%) 
IRAC galaxies assigned R^O.6 being the correct counter- 




Reliability 



Figure 6. The distribution of reliabilities for the IRAC candi- 
dates in the radio- and SMA-identified training set. The entire 
training set (black) is correctly identified (ranked first via the 
likelihood ratio; red) above i?=0.5, though there is contamination 
from possibly unassociated IRAC galaxies at almost all reliabil- 
ities (solid blue). The darker shaded region represents what we 
consider to be the robust range of reliabilities, while the lighter 
shaded region defines a tentative range. 

part. In GOODS-S, this ratio drops to 3/7 (43%) (see Ta- 
ble [TJ. 

4-2.2 Multiple Radio Systems and Radio Undetected 
SMGs 

The likelihood ratios and reliabilities for SMGs with mul- 
tiple radio counterparts can be seen in Table IA2I Some of 
these systems may be blended, multiple systems, and, as 
such, the IRAC counterparts for many of these systems have 
low reliabilities. There are a small number of cases (5/27), 
however, where a single IRAC candidate has an i?^0.8 and 
is thus likely associated with only one or none of the radio 
sources. One of these cases, AzTEC/C42, was observed with 
the SMA, which confirmed that the submm emission is only 
coming from one of the two radio sources within the AzTEC 
beam. The IRAC galaxy associated with the true radio coun- 
terpart was correctly identified with our technique. 

The ranking and reliabilities for the IRAC galaxies cor- 
responding to the radio-undetected SMG se10 can be seen 
in Table |M| There are 29 (70) radio-undetected SMGs with 
R^0.8 (R^0.6) (see Table [TJ. 



5 DISCUSSION 
5.1 Source Blending 

Multiple studies have established that some subset of SMGs 
are in fact blends of multiple sources of submm emission 

3 A list of AzTEC sources and potential IRAC 
counterparts, including photometry, is available at 
www. astro .umass . edu/aztec/CounterpartlD/counterpartid.html 
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which appear as one source in the large beams of single- 
dish submm observations. Our best estimates of the blending 
fraction come from interferometric studies, which have de- 
termined that 10-30% of known SMGs are mu l tiple systems 
I Younger et al.ll2007l . 120091 ; ISmolcic et ai1l2012l ; IBarger et all 
l2012h . A recent study of LABOCA sources imaged with 
ALMA determined that the overall blending fraction in their 
sample of 88 is 26%, though this fract ion is higher for t he 
very brightest SMGs (Saro^m >9 mJv: lKarim et al.ll2012h . 

The formulation of the reliability factor is such that it 
assumes there is only one unique counterpart. In the event 
of blending, the reliability of a potential counterpart will re- 
flect that there are multiple sources with a large LR values 
by assigning these sources correspondingly low reliabilities. 
These cases are indistinguishable from those in which an un- 
related IRAC galaxy is assigned a high LR due to a chance 
superposition and a redshift that is consistent with the peak 
in the IRAC color-color distribution of SMGs. For this rea- 
son, the reliability factor is biased against the identification 
of blended submm systems. 

Currently, there is no formulation of the reliability fac- 
tor that accounts for the possibility of multiple counter- 
p arts. An alternative to the reliability factor was presented 
in lChapin et al.l (|201ll ). which used simulations to determine 
a cut-off threshold on the likelihood ratio above which the 
rate of false positives was ~10 per cent. We performed a pre- 
liminary analysis to determine the utility of this technique 
with our dataset as follows: in COSMOS and GOODS-N, 
we choose 10,000 random locations with the areas around 
known SMGs masked out and calculate the LR of all IRAC 
galaxies within a search radius as described in Section [3] 
As the LR incorporates positional information based on the 
S/N of the submm observation, we repeated this process 
for a range of S/N bins. A likelihood ratio threshold is de- 
termined such that 90 per cent of the IRAC galaxies are 
rejected, which corresponds to a 10 per cent false positive 
rate. The LR of each potential IRAC counterpart is then di- 
vided by the appropriate thres hold according to the S/N of 
the corresponding SMC. As in lChapin et al.l (|201lT ). a nor- 
malized likelihood ratio (nLR) of greater than one indicates 
a candidate counterpart. 

An analysis of the results of using a nLR finds a higher 
fraction of SMGs with multiple sources assigned nLR>l, 
indicating blends, than we would expect from interfero- 
metric studies. In COSMOS, we find that 31/68 (46%) of 
SMGs have more than one counterpart with nLR>l and in 
GOODS-N this number is 8/13 (62%). In addition, we ex- 
amined m ultiple AzTEC source s that correspond to SCUBA 
sources m IBarger et all l|2012h . which were determined to 
not be blended at the resolution of the SMA. All five sources 
have multiple counterparts with nLR>l. 

These results indicate that the normalized likelihood 
ratio, when used in conjunction with our dataset and meth- 
ods, suffers from contamination from unrelated IRAC galax- 
ies. An additional difficulty with these simulations arises in 
small fields such as GOODS-N, where there are only ~100 
unique locations in the map given our typical search radius. 
While the normalized likelihood ratio is not biased against 
blended systems, it presents cases where blending is indis- 
tinguishable from confusion due to unrelated IRAC galaxies 
in excess of the 10% false positive rate expected. As such, we 
find the reliability factor to be more conservative in identi- 
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Figure 7. The dependence of the reliability factor on the AzTEC 
flux density for training set counterparts that are 1) ranked first 
(identified correctly) with our technique (filled symbols) and 2) 
those that are not ranked first (open symbols). There is no cor- 
relation, which suggests that the AzTEC flux density does not 
directly affect our counterpart identification technique. 



fying unambiguous, single counterparts and do not present 
the normalized likelihood ratios in this work. 



5.2 Dependence on the 1.1 mm and IRAC 
Observations 

In the previous section, we found that the accuracy of our 
technique, as blind-tested by our training set, improves with 
decreasing beamsize as we would expect. Here we explore ad- 
ditional effects on this technique; particularly we look at dif- 
ferences between the COSMOS and GOODS-S fields, since, 
despite being observed with similar beamsizes, our technique 
performs more accurately in COSMOS over GOODS-S by 
~10 per cent. 

The submm number counts reveal that SMGs in 
GOODS-S are sy stematically fainter than in COSMOS 
l|Scott et al.ll2012f ). Since our technique depends on the de- 
tectability of SMGs in all IRAC bands, we consider whether 
fainter submm sources implicitly translate to fainter IRAC 
counterparts. Our training set indicates that this is not 
the case as we find no correlation between the 1.1 mm 
flux density and any of the IRAC bands. This is what we 
would expect, given that the 1.1 mm flux density is pow- 
ered by dust emission and the IRAC bands primarily probe 
the mostly dust-free stellar bump at redshifts characteristic 
of SMGs. Additionally, we considered whether brighter or 
fainter SMGs have more easily identifiable counterparts due 
to some external factor which affects their IRAC colors, such 
as their redshift distribution. Figure [7] shows the reliabilities 
of SMGs in our training set versus their 1.1 mm flux den- 
sity; the two variables have a Pearson correlation coefficient 
of ~0.05, consistent with being uncorrelated. Cutting the 
GOODS-S AzTEC catalog to the depth of the COSMOS 
AzTEC catalog similarly produces no improvement in the 
counterpart identification results in GOODS-S. 

Next, we consider what effect the depth of the IRAC 
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catalog may have on our analysis. To directly compare 
GOODS-S and COSMOS, we cut the GOODS-S IRAC cat- 
alog to the depth of the COSMOS catalog, removing any 
sources with a 3.6 /jm detection less than 1 /xjy. This re- 
sults in an improved identification rate, with 7/15 (46%) 
IRAC counterparts to SMGs in GOODS-S correctly identi- 
fied. This is similar to the 49 per cent identification rate in 
COSMOS. This effect is due to the large beamsize of the 
AzTEC observations, necessitating a larger search radius 
for counterparts, and the depth of the IRAC observations 
in GOODS-S increasing the chances of detecting a high red- 
shift galaxy which is unrelated to the SMG, but, due to its 
redshift, has similar IRAC colors. These chance associations 
can be mitigated by higher resolution submm/mm obser- 
vations or by cutting the IRAC catalog, as SMGs will be 
brighter in the IRAC bands than other typical high redshift 
galaxies due to having larger stellar masses by about an or- 
der of magnitude (see lHainline et al.l l201ll . and references 
therein). 

5.3 Dependence on the Radio Flux Density 

The utility of our technique depends in large part on the 
applicability of the IRAC CDD to radio-undetected SMGs. 
We test this in two ways: first, we check that our IRAC color 
distribution is robust against the depth of the radio observa- 
tions. Second, we split our training set into bright and faint 
radio counterparts, to determine if there is a systematic dif- 
ference in their IRAC colors. 

The effects of the radio depth on our technique is tested 
by extending the COSM OS radio catalog from the origi- 
nal 4 . 5a cutoff to 3<r (e.g. IChapman et alJ r2005: Yu n et aT] 
2012). Since the false positive rate increases dramatically 
with decreasing S/N, we identify all ^3cr peaks above the 
local noise in the radio map and then match these peaks to 
the COSMOS IRAC catalog. Only those radio peaks with a 
corresponding IRAC counterpart within 2" are considered. 

Originally, we found 74 single radio counterparts to 
SMGs, plus 15 multiple radio systems. Using the deeper 
radio catalog, we find 26 new single radio detections and 
3 new multiple radio detections around previously radio- 
undetected SMGs. In addition, 7/74 COSMOS SMGs in 
the original training set were found to have a second ra- 
dio source within their search radius and one double sys- 
tem was found to have a third radio source (AzTEC/C184). 
This makes a total of 118/189 (62%) radio-detected SMGs in 
COSMOS, including 25 (13%) multiple radio systems (Ta- 
bles IAHA2|) . The IRAC colors of the new single radio de- 
tections can be seen overlaid on our original IRAC coun- 
terpa rt set in Figure [8] A 2-D K olmogorov-Smirnov (KS) 
Test (|Fasano &: Franceschinil ll987l returns a value of 0.84, 
indicating the original training set and extended radio coun- 
terparts are consistent with being drawn from the same par- 
ent distribution. When analysed using the IRAC CDD built 
with the original training set, 12/26 (46%) of the new radio 
sources were ranked first (4 with R^0.8). 

Combining the new single radio counterparts with the 
previous training set (and removing the new multiple radio 
systems), we perform the same kernel density estimation as 
before and apply this new IRAC color-color CDD in a blind 
test as in Section [4.11 We find that the addition of 26 new 
radio counterparts has a negligible effect on the resulting 
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Figure 8. IRAC colors of COSMOS radio sources in the extended 
catalog (S/N=3-4.5; red filled circles) overlaid on the original 
training set used in this study (see Figure [TJ. 



Table 2. Summary of Cryogenic Mis- 
sion Spitzer Legacy Programs 



Survey 


Area 


Reference 


GOODS 


280 arcmin 2 


Dickinson ("2003) 


SIMPLE 


-0.5 deg 2 


Damen et al. (20111 


SpUDS 


1 deg 2 


Dunloo et al. (20071 


S-COSMOS 


2 deg 2 


Sanders et al. (20071 


SDWFS 


10 deg 2 


Ashbv et al. (20091 


SWIRE 


~60 deg 2 


Lonsdale et al. (20031 



likelihood ratios and reliabilities assigned to each potential 
IRAC counterpart. This indicates that our IRAC color-color 
CDD is well-sampled and robust against the depth of the 
radio observations. 

As an additional check against potential bias on the 
IRAC colors used in this technique, we split our training 
set in half to analyse the IRAC colors of radio counterparts 
with radio flux densities greater than and less than 65 fj,3y. A 
2-D KS test of the two resulting distributions in IRAC color- 
color space gives a value of 0.98, confirming that the IRAC 
colors of SMGs with radio flux densities above and below 
65 /iJy are consistent with being drawn from the same parent 
distribution. This indicates that there is no dependence on 
the IRAC characteristic density distribution on the radio 
flux density. 

Given these two tests, we see no evidence that the IRAC 
colors of bright and faint radio counterparts are systemat- 
ically different, which is consistent with the detectability 
of radio counterparts being a function of the depth of the 
radio observations and a selection effect due to the strong 
positive k-correction in the radio regime. This result gives us 
confidence that the IRAC color-color CDD developed in this 
study can be applied to our radio-undetected SMG sample. 
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5.4 In the Era of Warm Spitzer, WISE, and JWST 

In this study, we have utilized data from the four IRAC 
bands available during Spitzer 1 s cryogenic mission. This 
mission surveyed 10s of square degrees on the sky over sev- 
eral wide, deep Legacy Programs, a summary of which can 
be seen in Table [2] These fields will be natural targets for 
future submm and mm observations. 

In May 2009, the cryogenic mission of Spitzer ended 
and it began the "warm" mission phase which includes only 
the 3.6 and 4.5 /im channels. Unfortunately, the close spac- 
ing of the two warm IRAC channels means that at z>2, both 
bands fall on the same side of the stellar bump and their 
redshift discriminating power is much diminished. As such, 
performing a characteristic density distribution analysis as 
seen here is not beneficial for SMGs discovered in fields with 
only warm Spitzer coverage, and we must resort to simpler 
selection methods, such as adopting a straight color cut or 
using a col or cut in the formation of an IRAC color p rior or 
P-statistic (|Yun et alJIiooi . 120121 : IChapin et al.ll201l| j. 

The Wide-fie ld Infrared Survey Explorer (WISE; 
IWright et ai]|2010h is an infrared survey which will map the 
entire sky at 3.4, 4.6, 12, and 22 /xm to 5a depths of 0.08, 
0.11, 1, and 6 mjy. Though this is far too shallow to detect 
the infrared counterparts of normal SMGs, WISE will be 
able to detect strongly lensed submm background galaxies , 
with magnification factors up to 30x (Eg ami et alJl2010tl. 



Thou g h these lensed source s are rare (e.g. Swinbank et al.l 



l20ld : iNeerello et ail l2010l: IConlev et a .1 l201lll. 



s urveys 



with Herschel (e .g. lEgami et alJ 120101: lEales et all l2010l : 
lOliver et all 120101), th e SPT (jVieira et alj|20ich . and ACT 
( Marriage et al.ll201ll ) and future large-scale SMG surveys 
with the Large Millimeter Telescope (LMT) and SCUBA- 
2 will dramatically increase the number of known lensed 
submm galaxies. The three short WISE bands will bracket 
the stellar bump at z>2 and have counterpart discriminat- 
ing power similar to that derived in this study. 

Hubble's successor, the James Webb Space Telescope 
(JWST), will have two instruments, NIRCam and MIRI, 
which together will provide imaging in the range 0.6 to 
28 (jm, covering the wavelength regime of Spitzer IRAC and 
MIPS 24 fim. With sensitivities ~50x better than IRAC, 
JWST will be able to detect mid-infrared counterparts out 
to very high redshifts in multiple bands which bracket the 
stellar bump. 



5.5 Use in Ongoing and Future Submm/mm 
Surveys 

We are currently seeing major advances in single-dish 
submm/mm observations, with improved resolution for 
observations at 450-1100 /jm. Large surveys are being 
carried out by SCU BA-2 on the JCMT (450, 850 /um; 
iDempsev et al.ll2012t ) and, in the near future, by AzTEC on 
the LMT (1100 /zm). To put our technique in the context 
of these surveys, we simulate the counterpart identification 
results that can be expected with the improved positional ac- 
curacy of these observations. We test two beamsizes: 14.5", 
which is the approximate FWHM of SCUBA-2 at 850 fxm, 
and 8.5", the FWHM of AzTEC on the LMT in its early 
32-m configuration. The latter beamsize is also close to the 
FWHM of SCUBA-2 at 450 /mi (~7"). Additionally, we sim- 
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Figure 9. The simulated percentage of correctly identified IRAC 
counterparts in COSMOS for beamsizes of 33" (circles), 14.5" (di- 
amonds) and 8.5" (squares), representing AzTEC on ASTE (this 
study), SCUBA-2 on JCMT, and AzTEC on LMT, respectively. 
The results using this technique (empty symbols) are an improve- 
ment over using just positions (tilled symbols) in the lowest S/N 
bins for all beamsizes. 



ulate 33", the beamsize of AzTEC on ASTE and the COS- 
MOS observations used in this study, which is comparable 
to the beamsize for SPIRE 500 /im observations. 

The simulations are done using the 68 SMGs in the 
COSMOS field that have a known radio and IRAC coun- 
terpart. To simulate the effects of random noise in the mm 
observations, a random position is drawn for each SMG from 
a distribution given by EqnQ] which is centered on the radio 
counterpart and has a FWHM determined by the observa- 
tions being simulated and a fixed S/N (Eqn[2]|. This random 
position determines the positional uncertainty, which is then 
combined with the IRAC colors of all candidate counterparts 
within an appropriate search radius as per the same blind- 
test performed in Section 14.11 This procedure is repeated 
1,000 times for each SMG for different S/N bins. The re- 
sults of the simulation can be seen in Figure [9j contrasted 
with the results of using positional information only to pre- 
dict the counterpart. Our technique is an improvement over 
using only position in the lowest S/N bins in all cases, and 
as expected the amount of improvement is largest for in- 
creasing beamsize. These simulations show that, even with 
higher resolution in future wide-area surveys, this technique 
can provide better predictions of SMG counterparts over po- 
sitional coincidence, as well as a measure of robustness of the 
prediction through the reliability factor. 



6 SUMMARY AND CONCLUSIONS 

In this study, we have developed a new technique for multi- 
wavelength counterpart identification to submm detections 
by expanding on previous color cut techniques to create a 
characteristic density distribution of IRAC colors for the 
SMG population. Our IRAC color-color distribution is based 
on a statistically significant set of 102 SMG counterparts 
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determined through radio and SMA detections and built 
using kernel density estimation with cross-validation. When 
combined with positional uncertainty information via the 
likelihood ratio, our technique provides a method for ranking 
IRAC galaxies surrounding SMGs according to their relative 
probability of being the correct counterpart. In addition, 
we use the reliability factor to weigh the importance of a 
given likelihood ratio against the likelihood ratios of other 
candidate IRAC counterparts. 

Once the IRAC color distribution is determined through 
known counterparts as in this study, it can then be applied 
to SMGs that lack other ancillary data, such as deep radio 
or MIPS 24 /im observations. The strength of the CDD for 
SMGs stems from the ubiquitous nature of IRAC detections 
of SMGs and the discriminating power of the IRAC colors, 
which can be used as a redshift indicator for galaxy popula- 
tions which 1) display a stellar bump in their infrared SEDs 
and 2) have a redshift distribution which peaks above z>l. 
Our overall technique, however, is general and can be applied 
to any population which occupies a unique color-magnitude 
or color-color space. 

In a blind-test using our training set, we show that our 
technique works well for mm observations with a moder- 
ate beamsize (e.g. 18" in GOODS-N or 14.5" for SCUBA-2 
surveys) with an identification rate of 85 per cent, compara- 
ble to the current identification rate in GO ODS-N, which 
utiliz es extremely deep radio observations (|Chapin et al.l 
2009). Our simulations demonstrate that this technique af- 
fords improvement over using only positional information, 
even for smaller beamsizes (e.g. 8.5" for AzTEC on the 
LMT). For larger beamsizes (>30"), increased positional un- 
certainties and chance associations with high redshift non- 
SMG galaxies reduce the identification rate to 33-49 per 
cent. These rates are comparable to current identification 
rates in COSM OS and GOODS-S l|Aretxaga et all l201ll ; 
lYunet alJlioTl ), using moderately deep radio observations. 
Our technique, however, identifies a different sub-population 
of counterparts than previous radio studies by using IRAC 
data alone to determine counterparts. These counterparts 
do not necessarily overlap with those identified through 
radio surveys and may have different properties, such as 
redshift distribution. Testing reveals no dependence of this 
method on the radio flux density, indicating it is applicable 
to radio-undetected SMGs. Future expansions of this tech- 
nique would include incorporating upper limits in order to 
allow for IRAC sources that are not detected in all four 
bands, and the expansion of the CDD beyond two dimen- 
sions to include additional information from a wider range 
of flux densities and colours. 

Though the challenges put forth by submm counter- 
part identification will be eased in the near future by in- 
struments and facilities such as SCUBA-2 on the JCMT, 
AzTEC on the LMT, and the relatively rapid interferome- 
try of ALMA, there will still be a great need for practical 
and efficient counterpart identification techniques. Single- 
dish facilities will increase the number of submm sources 
available for study by orders of magnitude within the next 
few years, making interferometric follow-up of each source 
impractical, while the need for multi-wavelength analy- 
sis remains vital for a complete picture of star forming 
galaxies and galaxy evolution. While this study only ad- 
dresses long- wavelength submm/mm observations, observa- 



tions at shorter- wavelengths in the FIR/submm such as with 
BLAST or Herschel suffer from similar counterpart identi- 
fication challenges, indicating the continuing importance of 
counterpart identification studies. 
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APPENDIX A: DATA TABLES 



Table Al: Training and Radio Counterpart Set: this table contains the radio counterparts to AzTEC SMGs in GOODS-N, 
GOODS-S, and COSMOS. The full table (385 rows) is available online. Sources are included in the training set if the radio 
counterpart was detected in the original radio catalogs (non-extended, see Section l5.3H and if there is a corresponding IRAC 
detection to the radio counterpart that is detected in all four bands. For COSMOS, radio detections which were found in the 
extended catalog (see Section 15.311 are indicated. Three additional radio- undetected SMA-confirmed IRAC counterparts are 
listed and included in the training set. For each source, the IRAC counterpart (if available) is listed on the same line as the 
radio counterpart, along with its rank, likelihood ratio, and reliability factor. The next two highest ranked IRAC galaxies 
are additionally listed (if available). Some AzTEC sources in COSMOS are cross-listed with Table lA2l due to having an 
additional nearby radio source in the extended radio catalog, making them multiple systems (see Section 15.311 . 



AzTEC ID a 


Rs 
(") 


Radio Coords b 
(J2000) 


P-stat 14GHz b 


IRAC Coords b 
(J2000) 


Rank b 


Likelihood 10 
Ratio 


Reliability 10 Other IDs c 
Factor 










GOODS-N 








AzTEC/GNl 


2.9 


J123711. 89+622211. 8 


0.002 


J123711. 88+622212. 5 

— 


1/1 

— 


182.47+4.7 

— 


1.00 


AzTEC/GN3 


5.2 


J123633. 42+621408. 7 


0.001 


J123633.40+621408.8 


1/4 


67.17+7.4 


.< 3 GOODS850-5 










J123633. 25+621411.4 


2/4 


20.49+5.8 


0.19 










J123633. 81+621407.5 


3/4 


16.40+10.3 


0.15 


AzTEC/GN4 


5.9 


J123550. 22+621041. 3 


0.014 


J123550. 35+621042. 2 


1/2 


152.65+3.8 


0.88 










J123549. 89+621043.1 


2/2 


21.22+3.8 


0.12 


AzTEC/GN5 


6.2 


J123730. 78+621258. 7 


0.010 


— 

J123730. 75+621258. 8 


— 
1/4 


— 

116.04+2.6 


— 

0.81 GOODS850-6 










J123730.08+621255.5 


2/4 


17.71+2.8 


0.12 










J123731. 09+621254.1 


3/4 


9.18+12.0 


0.06 


AzTEC/GN6 


6.3 


J123627.23+620605.7 


0.003 


J123627.21+620606.0 


1/5 


66.71+7.1 


0.65 










J123627.00+620610.0 


2/5 


15.92+18.7 


0.15 










J123627.20+620602.5 


3/5 


9.28+27.8 


0.09 


AzTEC/GN9 


7.2 


J123738. 15+621737.1 


0.004 


J123738. 26+621736. 8 


1/5 


49.92+4.6 


0.88 










J123737.54+621737.3 


2/5 


3.55+2.7 


0.06 










J123739.18+621736.5 


3/5 


2.49+0.8 


0.04 


AzTEC/GN12 


7.6 


J123632. 65+620621.1 


0.032 


J123632. 66+620621. 7 


1/2 


35.22+2.8 


0.95 










J123633.96+620613.8 


2/2 


1.81+3.4 


0.05 


AzTEC/GN14 


7.8 


J123651. 72+621221. 4 h 


0.031 


— 

J123651. 73+621220.7 


— 

2/5 


- 

8.54+0.5 


- 

0.39 GOODS850-1 










J123652. 08+621226.6 


1/5 


9.08+0.3 


0.41 










J123652. 69+621220.0 


3/5 


3.52+2.0 


0.16 


AzTEC/GN15 


7.8 


J123547.93+621529.2 


0.036 


- 


- 


- 


- 


AzTEC/GN19 


8.1 


J123604.40+620702.7 


0.010 










AzTEC/GN20 


8.1 


J123712. 48+621035. 4 


0.013 


J123712. 51+621036.0 


1/6 


141.25+4.0 


0.65 










J123712.49+621039.8 


2/6 


66.40+37.3 


0.31 










J12371 1.75+621036. 6 


3/6 


5.13+5.2 


0.02 


AzTEC/GN23 


8.6 


J123716. 67+621733. 2 


0.003 


J123716. 67+621733. 6 


1/8 


63.67+0.6 


0.51 










J123716. 14+621735.4 


2/8 


27.95+0.7 


0.22 










J123716. 16+621730.6 


3/8 


21.40+2.6 


0.17 


AzTEC/GN25 


8.7 


J123651. 72+620503.0 


0.024 










AzTEC/GN26 


8.7 


J123713. 86+621826. 2 


0.000 


J123713. 86+621826. 6 


3/7 


8.13+0.3 


0.09 










J123714.29+621824.3 


1/7 


51.98+19.1 


0.57 










J123713. 35+621830.3 


2/7 


17.04+23.1 


0.19 


AzTEC/GN29 


9.1 


J123619.89+621910.3 


0.109 










AzTEC/GN33 


9.2 


J123650.69+622016.5 


0.107 


J123650.68+622017.2 


2/8 


8.42+0.7 


0.23 










J123651. 01+622021. 7 


1/8 


12.90+14.2 


0.35 










J123650.38+622021.5 


3/8 


6.09+3.7 


0.17 


AzTEC/GN34 


9.3 


J123648. 21+622106. 9 


0.004 


J123648. 28+622107.2 


1/3 


120.88+3.6 


0.83 










J123647.37+622110.5 


2/3 


14.25+2.9 


0.10 










J123649. 42+622104.7 


3/3 


10.36+5.5 


0.07 


AzTEC/GN35 


9.3 


J123817.06+621426.3 


0.104 


















GOODS-S 








AzTEC/GSl 


5.0 


J033211. 37-275212.1 


0.026 


J033211. 35-275213.0 


1/1 


20.73+0.9 


0.99 
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Submm/mm Galaxy Counterpart Identification 17 



AzTEC ID a R s Radio Coords b P 1 4GHz b IRAC Coords b Rank b Likelihood Reliability Other IDs c 

(") (J2000) ' (J2000) Ratio b Factor b 



AzTEC/GS5 


7.5 


J033151. 11-274437. 6 


0.048 


J033151. 09-274437.0 


2/3 


3.04+0.1 


0.33 










J033151. 06-274432. 4 


1/3 


5 58+4 

7 a '.V 7/ 1_ i a \ ) 


0.60 










J033151. 47-274427. 9 


3/3 


0.41+0.1 


0.04 


AzTEC /GS7 


8.5 


1033213 85-275600 3 


0.079 


103321 3 85-275559 9 

U \JyJyj£j±*-J*yJyJ ^il JUJi/iC 


4/5 


4.10+0.2 


0.03 










J033213. 56-275608. 3 


1/5 


100.79+37.5 


0.83 










1033213 91-275607 4 


2/5 


8 22+3 6 


0.07 


AzTEC/GS8 


8.6 


J033204. 90-274647. 5 


0.023 


J033204. 87-274647. 3 


2/4 


26.61+0.8 


0.34 










1033205 34-274643 9 


1 / 1 
X/<± 


48.54+2.2 


0.62 










J033205. 29-274639.1 


3/4 


2.73+0.1 


0.03 


AzTEC/GS9 


8.5 


J033303. 02-275 146. 5 


0.043 


J033303. 00-275146. 3 


1/5 

X/ *J 


31 69+2 


0.41 










1033302 91-275151 


2/5 


oc; 24+6 6 


0.32 










1033302 45-275145 3 


3/5 


20 31+7 8 


0.26 


AzTEC /GS10 


8.9 


7033207 31-275120 8 


0.034 


J033207. 28-275120. 2 


2/3 


7 26+0 3 


0.18 










1033207 09-275129 


1/3 


31 71 + 14 


0.79 










1033207 34-275134 4 

U V/'7<7^jV7 . <7 i — / 7 1 './i . i 


3/3 


0.94+1.4 


0.02 


AzTEC/GS20 


10.3 


J033235. 08-275532. 6 


0.025 


J033235. 07-275532. 6 


1/3 

X/ o 


8.87+0.0 


0.72 










1033234 63-275542 8 


2/3 


2 39+0 2 

^— p ' } 1_ V 7 p £t 


0.19 










1033235 60-275543 7 


3/3 


92+0 2 

V7 . v 7 — 1_ V7 . »£i 


0.07 


AzTEC/GS21 


11.0 


T033247 59-274452 3 


0.062 


J033247. 55-274452. 4 


3/6 


9.84+0.4 


0.18 










1033247 94-274450 1 


1/6 


28 87+10 5 


0.53 










J033247. 24-274444. 4 


2/6 


1 1 39+0 9 


0.21 


AzTEC/GS25 


13.0 


J033246. 85-275121.0 


0.053 


1033246 78-275120 9 


1/6 


25.23+0.8 


0.85 










1033246 48-275128 2 


2/6 


q ns+2 i 


0.10 










J033247. 67-275129.0 


3/6 


1.09+4.8 


0.04 


AzTEC /CS31 


13.4 


1033243 48-273929 4 


0.069 


1033243 49-273929 1 


3/7 
°/ ' 


Zj . -Mr 1 I 1 


0.22 










J033243. 32-273924.0 


1/7 


3.29+22.4 


0.28 










1033242 80-273927 1 


2/7 


2.94+0.1 


0.25 


AzTEC /GS33 


12.7 


T033248 79-275314 8 


0.061 


1033248 77-275314 4 


2/8 


26 71 +4 3 


0.36 










J033249. 71-27531 1.2 


1/8 
x/o 


41 07+16 4 


0.56 










1033249 58-275322 6 

•J V ' '7 <7 £t i \ ' . 7 >7 — f • 7 ■ / — Z_ . V7 


3/8 


3 58+3 1 

*-J p f7 O 1_ t7 p -L 


0.05 


AzTEC /CS34 


13.9 


1033229 96-274301 8 


0.140 


1033229 93-274301 


9 /12 




0.02 










J033230. 06-274306.7 


1/12 


7.83+1.2 


0.28 










J033229. 68-274306. 3 


2/12 


6.22+0.5 


0.22 


AzTEC/GS38 


14.2 


J033209.72-274248.6 


0.069 


J033209. 71-274248. 3 


8/8 


1.21+0.1 


0.03 










J033208. 74-274248.6 


1/8 


12.30+0.5 


0.32 










J033209. 69-274238. 7 


2/8 


7.59+2.6 


0.20 


AzTEC/GS51 


15.3 


J033157.59-274503.8 


0.325 


J033157.60-274504.1 


2/7 


22.52+2.1 


0.27 










J033157. 72-274510.1 


1/7 


28.78+15.5 


0.35 










J033156. 97-274501. 2 


3/7 


15.32+8.6 


0.19 


AzTEC/GS54 


16.0 


J033243. 20-275514. 3 


0.004 


J033243. 19-275514.4 


1/12 


12.94+0.3 


0.32 










J033243. 75-275516. 2 


2/12 


11.31+1.5 


0.28 










J033243. 23-275509.2 


3/12 


3.01+12.5 


0.08 










COSMOS 








AzTEC/Cl d 


6.6 


J 100141. 75+022713. 4 


0.002 


J 100141. 86+022713. 7 


1/2 

X/ ^ 


8.82+0.7 


0.75 










J 100141. 51+022714. 7 


2/2 


2.71+0.2 


0.23 


AzTEC/C2 


6.2 














AzTEC/C3 f 


6.6 


J 100008. 02+02261 2.1 


0.003 










AzTEC/C4 g 


6.6 






J095931. 69+023043. 9 


1/2 


68.98+8.7 


0.86 










J095931. 99+023042. 5 


2/2 


11.14+5.2 


0.14 


AzTEC/C6 


7.1 


J100057.26+022012.2 


0.007 










AzTEC/C8 d 


7.1 


J100013. 82+015639. 5 


0.023 


J100013. 85+015639.4 


2/3 


11.02+0.9 


0.27 










J100014.30+015641.5 


1/3 


27.01+7.4 


0.67 










J100014.63+015632.5 


3/3 


2.14+0.9 


0.05 


AzTEC/C9 


8.2 


J100122. 97+022005.4 


0.012 


J100122. 95+022005. 9 


2/5 


42.69+3.9 


0.31 










J100123.15+022004.9 


1/5 


79.42+8.1 


0.58 










J100123.43+022001.6 


3/5 


10.52+6.0 


0.08 


AzTEC/C10 f 


8.3 














AzTEC/C12 d 


8.2 


J100136. 78+021110. 2 


0.026 


J100136. 81+021110.0 


2/2 


14.38+2.5 


0.35 



AzTEC8 
AzTEC2 
AzTECM 



AzTEC15 
COSLA-17 
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18 Alberts et al. 



AzTEC ID a 


Rs 


Radio Coords b 


P i.4GHz b 


IRAC Coords b 


Rank b 


Likelihood 


Reliability 


Other IDs c 




(") 


(J2000) 


(J2000) 




Ratio b 


Factor b 












J100136. 61+021058.7 


1/2 


25.87+7.1 


0.64 




AzTEC/C13 


9.1 


J095837.96+021408.5 


0.000 


J095837.96+021408.4 


1/7 


61.08+3.1 


0.39 


COSLA-54 










J095837.76+021402.7 


2/7 


46.54+6.8 


0.30 












J095837.80+021412.5 


3/7 


36.09+15.2 


0.23 




AzTEC/C14 


9.8 


J095957.29+022730.4 


0.002 


J095957.29+022730.6 f 








AzTEC9 


AzTEC/C15 


10.0 


J100132. 00+022515.0 


0.030 


J100132. 03+022514. 7 


2/5 


3.29+0.1 


0.24 












J100131. 50+022516.1 


1/5 


5.46+4.6 


0.40 












J100132.03+022514.7 


2/5 


3.29+0.1 


0.24 




AzTEC/C16 


10.0 


J095854.17+021646.0 


0.010 


J095854. 20+021646.0 


1/6 


67.38+5.8 


0.61 












J095853. 69+021652.9 


2/6 


33.11+3.4 


0.30 












J095854.37+021653.5 


3/6 


5.31+1.1 


0.05 




AzTEC/C17 


10.4 


J100055. 37+023441. 5 


0.040 


J100055. 39+023441. 5 


1/2 


8.75+0.1 


0.61 












J100055.03+023435.7 


2/2 


5.34+0.6 


0.37 




AzTEC/C18 


10.6 


J100035. 29+024353. 3 


0.006 


J100035. 30+024353. 3 


1/3 


100.09+3.4 


0.74 


AzTEC12 










J100034.98+024349.5 


2/3 


28.94+7.6 


0.21 












J100034.83+024400.3 


3/3 


6.09+2.4 


0.04 




AzTEC/C20 d 


9.1 


J100114.54+022705.4 


0.005 


J100114. 54+022705.4 


1/5 


91.81+14.6 


0.79 












J100114. 78+022703.6 


2/5 


20.25+10.2 


0.17 












J100114.28+022655.4 


3/5 


2.28+0.4 


0.02 




AzTEC/C22 


11.0 


J100008. 91+024010. 3 


0.008 


J100008. 94+024010. 6 


1/4 


59.51+3.6 


0.95 


AzTEC 11 










J100008.56+024013.3 


2/4 


1.47+0.9 


0.02 












J100008.75+024018.7 


3/4 


1.12+0.3 


0.02 




AzTEC/C24 


11.0 


J100010. 36+022224.8 


0.008 


J100010. 36+022224. 5 


1/3 


52.86+5.4 


0.70 












J100010.08+022217.9 


2/3 


19.77+18.1 


0.26 












J100010. 70+022224.8 


3/3 


2.20+1.2 


0.03 




AzTEC/C25 


11.1 


J100121. 95+015644.1 


0.050 


J100121. 95+015644.1 


1/7 


8.54+0.3 


0.29 












J100121. 60+015646.1 


2/7 


7.66+1.7 


0.26 












J100121. 47+015656.6 


3/7 


6.75+2.3 


0.23 




AzTEC/C28 


11.3 


J095849. 28+021301. 8 


0.034 


J095849. 29+021301. 3 


2/7 


26.01+1.7 


0.33 












J095849. 27+021255.1 


1/7 


45.95+7.9 


0.57 












J095849.29+021301.3 


2/7 


26.01+1.7 


0.33 




AzTEC/C29 


11.3 


J095918. 37+020107.1 


0.001 


J095918. 37+020107.5 


2/3 


14.08+5.4 


0.17 












J095918.60+020104.5 


1/3 


55.39+30.9 


0.69 












J095918.37+020107.5 


2/3 


14.08+5.4 


0.17 




AzTEC/C31 


11.6 


J100147.01+022447.7 


0.050 


J100146. 97+022446. 9 


3/4 


2.66+0.2 


0.11 












J100147.72+022456.4 


1/4 


14.66+18.5 


0.63 












J100147. 19+022452.2 


2/4 


5.42+0.5 


0.23 




AzTEC/C32 


11.8 


J100012. 53+020124.5 


0.020 












AzTEC/C33 


11.8 


J100027.14+023140.6 


0.043 


J100027.15+023140.9 


2/6 


11.53+0.6 


0.32 












J100026.66+023127.2 


1/6 


13.10+20.8 


0.36 












J100027. 15+023140.9 


2/6 


11.53+0.6 


0.32 




AzTEC/C35 


12.0 


J100008. 97+022026. 5 


0.041 


J100008. 94+022026.1 


1/6 


30.73+2.9 


0.31 












J100008. 65+022030.5 


2/6 


27.20+12.4 


0.28 












J100007.92+022023.1 


3/6 


21.02+5.3 


0.21 




AzTEC/C36 


12.0 


J095840. 27+020514.7 


0.000 


J095840. 28+020514. 6 


3/3 


1.75+0.1 


0.04 












J095840.33+020519.7 


1/3 


43.04+18.7 


0.87 












J095839.79+020512.5 


2/3 


4.23+7.1 


0.09 




AzTEC/C37 


12.0 


J100121. 82+023129. 3 


0.031 


J100121. 87+023129. 3 


1/2 


18.89+4.7 


0.98 












J100121. 39+023123.4 


2/2 


0.02+13.2 


0.00 




AzTEC/C38 


12.0 


J100023. 67+022155. 3 


0.025 


J100023. 65+022155. 3 


1/3 


36.93+2.6 


0.93 


COSLA-38 










J100023.78+022148.7 


2/3 


2.16+4.8 


0.05 












J100023.24+022139.3 


3/3 


0.34+2.7 


0.01 




AzTEC/C40 


12.2 


J095935.34+021920.1 


0.063 


J095935. 36+021919.3 


2/3 


1.02+0.3 


0.15 












J095934.76+021930.2 


1/3 


4.78+3.4 


0.71 












J095935. 36+021919.3 


2/3 


1.02+0.3 


0.15 




AzTEC/C41 d 


9.8 


J100148. 19+022132. 8 


0.008 


J100148. 17+022132. 3 


1/4 


49.91+1.7 


0.92 












J100147.99+022129.1 


2/4 


2.92+0.5 


0.05 












J100147.30+022130.3 


3/4 


1.26+0.5 


0.02 




AzTEC/C42S 


12.5 


J100019. 74+023204.1 


0.00 


J100021. 13+013852. 6 


1/9 


21.10+2.3 


0.58 


AzTEC5 






J100019.98+023209.9 


0.028 


J100108. 35+015712.9 


2/9 


6.33+0.4 


0.17 
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Submm/mm Galaxy Counterpart Identification 19 



AzTEC ID a R s Radio Coords b P 1 4GHz b IRAC Coords b Rank b Likelihood Reliability Other IDs c 

(") (J2000) ■ - (J2000) Ratio b Factor b 











1100019 49+023211 


3/9 


3.18+1.7 


0.09 


AzTEC /C43 e 


12.5 


J 100003. 12+020201. 5 


0.032 


J 100003. 09+020201. 5 


5/8 


7.91±2.3 


0.05 










J 100003. 40+020204. 4 


1/8 


65.08±7.8 


0.38 










1100003 84+020207 2 


2/8 


47 53+12 1 


0.28 


AzTEC /C44 d 


10.0 


1100034 23+014858 2 


0.023 


J 100034. 25+014857. 8 


1/9 


38 60+8 9 


0.47 










J100033. 77+014852. 5 


2/9 


22.54±6.6 


0.27 










1100033 89+014911 

U \-\J\J\J\J*-J * UJ | u _L TEC -L _L » V7 


3/9 


7.81±0.2 


0.09 


AzTEC/C45 


12.6 


1100006 83+023307 9 


0.065 










AzTEC /C46 


12.6 


J 1001 14. 71+023518. 3 


0.021 


J100114. 74+0235 17. 5 


1/6 


87.64±14.0 


0.85 










1100114 17+023501 9 


2 t ) 


6 95+2 3 

V » . v 7 } 1_ — . • } 


0.07 










1100114 05+023504 7 

U -L\7V7-l-_l_Tt.l7(7 | U 'J U U T: ■ 1 


3/6 


4.07±1.6 


0.04 


AzTEC/C47 


12.6 


J095940. 87+0201 13. 3 


0.038 


J095940. 86+0201 13.3 


1/6 


9.71±0.7 


0.40 










1095941 56+020104 7 


2/6 


7.49±4.2 


0.31 










1095940 66+020103 7 


3/6 


K Z.A4-0 n 

'J . (Jri_l_.ii ."J 


0.23 


AzTEC /C48 


12.6 


J 100039. 28+023845. 3 


0.003 


J 100039. 29+023845. 2 


2/4 


40.31±3.1 


0.35 










1100039 78+023850 7 

U 1-\J\J\Jk)<j . 1 O |^ VJ Zj >7 O (7 W . 1 


1 /4 


41 04+1 2 

rt _L ■ _*rr_l L _. . 17 


0.35 










1100039 29+023845 2 

tJ 1 \ J \ 1 \ 1 ' J , J . — , / | 17 v7 *-7 (7 • _* 


2/4 


40 31+3 1 


0.35 


AzTEC/C50 


12.7 


J095933. 20+020832. 8 


0.001 


J095933. 27+020832. 7 


1/4 


61.21±9.0 


0.97 










1095933 61 +020836 9 


2/4 


1 00+2 

_L * \j \ J _1_ _j a VJ 


0.02 










1095932 75+020824 7 


3/4 


0.42±1.6 


0.01 


AzTEC/C52 


12.8 


J 100156. 57+022 101.1 


0.026 


J100156. 57+022100. 9 


2/3 


6.75±0.2 


0.22 










11 001 56 47+0221 07 4 


1 /3 


on 27+2 2 


0.65 










1100156 57+022100 9 


2/3 


6 75+0 2 


0.22 


AzTEC/C55 


12.9 


J100004. 41+015515. 3 


0.074 


J100004. 37+015516.1 


5/7 


3.50±0.3 


0.05 










1100005 09+015518 2 


1/7 


39 30+3 9 


0.58 










11 00004 82+01 5513 9 

U -L (7l7>_H7Tt. Oil 7^17 -L '.J'J -LrJ. 


2/7 


10.14±7.7 


0.15 


AzTEC/C57 


13.0 


J095958. 23+021314. 3 


0.049 


J095958. 23+02 1314. 3 


3/4 


0.80±0.1 


0.15 










1095958 49+021316 1 


1/4 


3 09+1 3 

ij p \J C _l L p 


0.57 










1095959 03+021323 3 


2/4 


80+0 1 

v7 * O \J 1_ \J p -L 


0.15 


AzTEC/C58 


13.2 


J 100020. 07+014502. 6 


0.019 


J100020. 06+014502. 9 


1/4 


43.44±6.3 


0.83 










1100020 86+014459 7 

U 1 V7 V7 \7 — V7 . \J VJ | V7 1 T_ r ' V - 7 . 1 


2/4 


5 86+2 

7 . v v_' i_ — p \ j 


0.11 










1100020 23+014510 6 

O _1_\7 W\7 — J V7 . — J i7 7^ Ul^'JlU.U 


3/4 


2.10±0.7 


0.04 


AzTEC/C61 


13.2 


J 100120. 06+023443. 7 


0.018 


J100120. 09+023443. 6 


2/6 


25.03±1.7 


0.33 










11 001 1 9 74+023434 7 


1 /fi 


Q2 69+1 6 5 


0.43 










1100120 09+023443 6 

U 1 ' ' V ' 1 _j V7 . V7 .7 | V7 — ' 1 r < 7 ■ vJ 


2/6 


25 03+1 7 


0.33 


AzTEC/C65 


13.4 


7095942 94+022145 


0.046 


J095942. 96+022145.0 


4/10 


4 36+0 2 

1 . ■ J v_' i_ v. / p — • 


0.14 










1095942 72+0221 38 


1/10 


g 34-I-4 7 


0.21 










1095942 76+022143 3 


2/10 


4 fi3+n 5 

t:.UO_LU. <-J 


0.15 


AzTEC/C66 


13.4 


J 100104. 64+022634.0 


0.015 


J 100104. 63+022633. 4 


1/4 


20.19±1.4 


0.60 










1100104 34+022639 5 


2 /4 


8.21±2.9 


0.24 










J100104. 57+022627. 2 


3/4 


3.47±0.5 


0.10 


AzTEC/C71 


13.6 


J095953. 86+021853. 8 


0.017 


J095953. 91+02 1853. 6 


1/5 

X/ o 


13.44-1-0.9 


0.51 










1095953 57+021854 3 


2/5 


S 49+0 8 


0.32 










1095954 19+021851 1 


3/5 


•3 01 +3 r> 


0.15 


AzTEC/C72 


13.6 


J100158. 98+020457. 8 


0.075 


J100159. 00+020457. 8 


2/3 


3.89±0.5 


0.23 










1100160 00+020509 4 

O 1UU1 UU . T^lJii'J'JiJC'.Tt 


1 /3 


1 2 fifi+fi 1 


0.74 










1100159 00+020457 8 


2/3 


3 89+0 5 

iJtUC _1_ \J p (7 


0.23 


AzTEC /C74 


13.8 


J 100105. 86+022 133.1 


0.025 


J 100105. 87+022132. 8 


3/5 


4.07±0.4 


0.14 










11 001 05 80+0221 36 

O 1 V7 V 1 \ > • 1 . 7" V7 | V7 — 1 ' ? V f . V ) 


1/5 

-7 u 


19 64+7 9 

_LC.V7^r_l_ 1 • i7 


0.69 










J100105. 33+022132. 5 


2/5 


4.17±4.5 


0.15 


AzTEC/C77 


13.9 


J095935.31+015759.3 


0.014 


J095935. 22+01 5759. 6 


1/4 

x/ 


20.68±6.4 


0.59 










1095934 89+015748 5 


2/4 


6 34+3 5 

V7 . ' J 1 1_ ' ' p (7 


0.18 










1095935 42+015754 


3/4 


4.11+1.2 


0.12 


AzTEC/C79 d 


10.0 


J095943. 62+021340. 4 


0.032 


J095943. 65+021340.1 


1/4 


37.98±2.1 


0.90 










J095944.31+021342.0 


2/4 


3.46±1.6 


0.08 










J095944.46+021355.0 


3/4 


0.30±0.1 


0.01 


AzTEC/C82 d 


10.4 


J100116. 68+021647.7 


0.003 


J100116. 63+021647.2 


2/12 


3.49±0.8 


0.15 










J100117.19+021645.5 


1/12 


3.85±0.3 


0.16 










J100117.45+021646.4 


3/12 


3.18±0.3 


0.14 


AzTEC/C84 


14.2 


J095942. 59+015501. 6 


0.051 


J095942. 57+015501.6 


2/3 


6.42±0.2 


0.17 










J095943. 38+015510.7 


1/3 


29.80±10.8 


0.79 










J095942. 57+015501.6 


2/3 


6.42±0.2 


0.17 


AzTEC/C85 


14.2 


J100139. 75+022548. 8 


0.039 


J100139. 76+022548. 5 


3/6 


1.20±0.0 


0.06 



© 2013 RAS, MNRAS 000.[TJl30l 



20 Alberts et al. 



AzTEC ID a R s 
(") 



Radio Coords b 
(J2000) 



P i.4GHz k 



IRAC Coords b 
(J2000) 



Rank b 



Likelihood 
Ratio b 



Reliability 
Factor b 



Other IDs c 



AzTEC/C86 14.2 J100109. 28+021721. 7 0.016 



AzTEC/C87 d 10.6 J100204.95+021701.9 0.035 



AzTEC/C88 d 11.0 J095937.44+020424.6 0.001 



AzTEC/C91 e 14.5 J100128. 20+022340. 8 0.038 



AzTEC/C93 14.6 J100131. 88+021138. 8 0.004 



AzTEC/C95 d 11.0 J100018. 19+021243.0 0.004 



AzTEC/C97 14.7 J100214.51+021942.9 0.013 



AzTEC/C98 14.7 J100043. 18+020519.0 0.005 



AzTEC/C103 14.8 J100123. 86+015613.1 0.036 



AzTEC/C105 d 11.1 J095845. 11+021431.3 0.058 



AzTEC/C109 14.9 J100111. 56+022840. 8 0.004 



AzTEC/Clll 15.0 J095929. 22+021243. 8 0.019 



AzTEC/C112 15.1 J100011. 04+015314.2 0.014 



AzTEC/C113 15.1 J095914.40+023000.1 0.036 



AzTEC/C115 d 11.3 J100015. 31+020532. 2 0.026 



AzTEC/C116 15.2 J100109. 86+020346. 4 0.008 



AzTEC/C118 15.2 J095959. 95+020633. 3 0.023 



AzTEC/C121 d 11.3 J095952. 08+020930. 5 0.020 



AzTEC/C123 d 11.6 J100022. 81+015140. 9 0.014 



AzTEC/C126 15.4 J100158. 81+022240. 2 0.048 



AzTEC/C127 15.5 J100125. 31+023527.5 0.007 



1100140 49+022542 fi 


1 /6 


io 35+0 3 


0.59 


1100139 76+022538 6 


2/6 


4.82+0.2 


0.23 


J 100109. 29+02 1721. 9 


1/5 


35.84+5.9 


0.82 


11 001 09 1 1 +021 720 


2/5 


4.05+6.4 


0.09 


11 001 09 91 +021 727 6 


3/5 


3 32+0 1 


0.08 


J 100204. 93+02 1701.0 


2/10 


15.83+1.0 


0.21 


1100205 48+021657 9 

U 1 V7 V7 V / ' 7 . i v 7 | V7 — 1 V7<7 1 . v ' 


1/10 


27 95+4 

J— 1 • • ./ " J 1_ i . V7 


0.38 


1100205 76+021659 4 


3/10 


10 73+0 4 

-L17 . I ' 7 1_ ^.7 • 


0.14 


J095937. 38+020424. 3 


3/6 


2.49+0.7 


0.19 


1095937 45+020434 3 

U V7 ,7 \J v ' ' ./ / . i ' ) \ V7 — \7 i '7 i ■ i } 


1/6 

1/ u 


3.13+4.7 


0.24 


1095937 38+020424 3 

U V./ >J "7 v' ' 7 / . <7 >7 | V7 — \7 i i . ' } 


3/6 


2.49+0.7 


0.19 


J 1001 28. 2 1+022340. 8 


3/8 


2.22+0.1 


0.07 


J100128. 48+022345.1 


1/8 


19.11+6.7 


0.57 


1100129 30+022348 2 


2/8 


1 • ' J !_» _l_ \J • <-J 


0.22 


J 100131. 87+02 1138. 7 


1/5 


31.70+1.5 


0.79 


11 001 32 04+021 1 44 9 


2/5 


4.02+6.1 


0.10 


11 001 32 58+021 1 30 7 


3/5 


2.07+2.7 


0.05 


J100018. 22+02 1242. 4 


2/6 


18 96+0 9 

ID* i7 V7 1_ V 7 . .7 


0.35 


J100018. 54+021241.1 


1/6 

1/ u 


29 06+18 4 


0.54 


1100017 82+021231 5 


3/6 


3.70+1.4 


0.07 


J100214. 49+021942. 7 


2/5 


32 09+1 8 

'7 *JiU7 _l L . O 


0.39 


1100214 71+021941 2 


1 /5 


47 nfi+S 9 

L ± 1 .1111 1 ( 1 ■ .~7 


0.57 


J100214. 49+021942. 7 


2/5 


32 09+1 8 


0.39 


J100043. 18+020519. 2 


1/6 


56.60+3.1 


0.47 


1100043 03+020516 8 


2/6 


S3 74+2 5 

■ ) '7 . f i ^ . '7 


0.28 


J100042. 48+020520. 8 


3/6 


16 33+2 9 


0.14 


J100123. 86+015613. 5 


6/8 


3.67+0.4 


0.03 


11 001 24 58+01 5609 8 


1 /8 




0.47 


11 001 24 48+01 561 6 

O i V7 *7 liJT.^U | V7 i '7 V7 -L vJ . W 


2/8 


33 61 +1 3 


0.29 


J095845. 13+021430.9 


1/7 


10.54+0.7 


0.41 


1095844 56+021443 6 


2/7 


3 95+0 4 


0.15 


1095845 33+021456 

•J V7 v7 "7 (7 i ' 7 . ' 7 ' } | V7 — 1 i ' 1 V7 • \7 


3/7 


3 91 +0 6 


0.15 


J10011 1.61+022841.0 


4/9 


38 03+2 4 

• J 1 - ' . V7 V7 1_ — • I 


0.17 


1100111 68+022832 7 


i /q 

x/a 


«c S6+1 5 2 


0.39 


J10011 1.33+022831. 8 


2/9 


47.21+9.9 


0.21 


T095929 21+021244 


1 /6 


66 1 1 +7 

UUpIJ 1 1 ■ • J 


0.87 


1095929 83+021247 7 


2/6 


3 52+2 

' 7 i " 7 i— . V7 


0.05 


1095929 47+021242 9 


3/6 


2.82+2.8 


0.04 


J10001 1.03+015314. 3 


1/5 

X/ u 


20.56+1.5 


0.82 


1100011 18+015258 3 


2/5 




0.11 


1100011 53+015311 8 

U X \7 17 V7 J- -1- . '7 * 7 | V/ J. '.J \J .L J. . C7 


3/5 


1.32+0.7 


0.05 


J095914. 42+022959. 8 


1/3 

X/ o 


21.70+2.0 


0.49 


1095915 25+022950 9 


2/3 


19 36+12 9 


0.44 


109591 5 60+023000 8 

JUi/JC'liJ. 17VJ |^ V7i-J '7V7V7 W . O 


3/3 


o SQ+n 3 


0.07 


J100015. 36+020531. 8 


1/4 

X/ 


37.58+8.6 


0.87 


1100014 80+020526 7 


2/4 


q 23+2 7 


0.07 


J100015. 59+020533. 9 


3/4 


1.97+3.1 


0.05 


J 100109. 83+020345. 9 


1/7 


73.97+12.9 


0.47 


1100110 09+020345 7 


2/7 


fi7 94+25 7 


0.43 


1100109 40+020353 1 

•7 1 17 V' 1 \ ' .7 . i V7 | V7 — V7 ' 7 ' ' v > . 1 


3/7 


9 86+10 5 


0.06 


J095959. 98+020634.0 


3/9 


15.26+1.8 


0.24 


109^959 69+020627 1 


1 /q 


i a fic;+i s n 


0.26 


J095959. 75+020634. 5 


2/9 


16.07+6.5 


0.25 


J095952. 09+020931. 9 


3/3 


0.05+0.0 


0.06 


JOyoyol. 6o+0209l0. / 


1 /*j 
1/3 


0.o/±0.1 


0.38 


J095952.09+020931.9 


3/3 


0.05+0.0 


0.06 


J100022. 82+015140. 6 


1/9 


30.72+2.6 


0.29 


J100023. 00+015144.0 


2/9 


27.87+2.5 


0.26 


J100022. 18+015152.8 


3/9 


21.24+8.0 


0.20 


J100158. 83+022240.1 


5/5 


1.95+0.1 


0.07 


J100159.01+022235.3 


1/5 


18.28+7.7 


0.61 


J100159.63+022236.0 


2/5 


3.52+0.9 


0.12 


J100125. 32+023527.2 


1/4 


59.01+2.3 


0.82 
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COSLA-18 



Submm/mm Galaxy Counterpart Identification 21 



AzTEC ID a 


Rs 


Radio Coords b 


P 1 ,PH7 b 

1 .4V_jrlZ 


IRAC Coords b 


Rank b 


Likelihood 


Rcliabil 




(") 


(J2000) 


(J2000) 




Ratio b 


Factoi 










J100125. 54+023530.6 


2/4 


4.91+0.9 


0.07 










J100125.70+023527.0 


3/4 


4.34+0.8 


0.06 


AzTEC/C130 


15.6 


J095957.35+020311.3 


0.018 


J095957.34+020311.6 


1/6 


33.23+2.3 


0.51 










J095956.82+020310.9 


2/6 


14.12+16.1 


0.22 










J095957.39+020259.5 


3/6 


13.73+2.8 


0.21 


AzTEC/C131 


15.7 


J100010.08+023745.7 


0.107 


J100010. 15+023745.3 


4/6 


4.14+0.1 


0.13 










J100011. 08+023749.4 


1/6 


11.59+1.0 


0.36 










J100011. 29+023801. 2 


2/6 


11.10+1.1 


0.34 


AzTEC/C133 


15.8 


J100120. 854-022623. 3 


0.006 


J100120. 84+022623. 5 


1/6 


48.28+5.1 


0.61 










J100119.97+022629.2 


2/6 


14.26+1.7 


0.18 










J100120.65+022623.6 


3/6 


11.03+6.3 


0.14 


AzTEC/C136 


15.9 


J095933. 474-022347.4 


0.001 


J095933.42+022346.8 


1/7 


43.51+3.2 


0.84 










J095933.84+022349.3 


2/7 


4.68+0.4 


0.09 










J095933. 06+022338.8 


3/7 


1.07+0.3 


0.02 


AzTEC/C137 d 


11.8 


J095953. 69+022233. 8 


0.040 


J095953. 71+022232. 9 


3/5 


3.23+0.2 


0.11 










J095953.28+022243.3 


1/5 


14.04+10.9 


0.49 










J095953.71+022232.9 


3/5 


3.23+0.2 


0.11 


AzTEC/C138 d 


11.8 


J100019.56+023511.3 


0.107 


J100019. 52+023511. 9 


3/5 


0.96+0.3 


0.13 










J100020.69+023520.5 


1/5 


2.81+3.6 


0.38 










J100019.52+023511.9 


3/5 


0.96+0.3 


0.13 


AzTEC/C140 c 


16.0 


J100124.96+015147.3 


0.004 


J100125. 00+015147.2 


2/5 


37.98+4.8 


0.35 










J100124.60+015141.3 


1/5 


61.43+9.8 


0.57 










J100125.00+015147.2 


2/5 


37.98+4.8 


0.35 


AzTEC/C142 e 


16.0 


J100018.82+020245.0 


0.070 


J100018.83+020245.0 


6/7 


0.93+0.0 


0.03 










J100018. 37+020245.1 


1/7 


18.97+14.7 


0.61 










J100017.47+020252.3 


2/7 


4.27+0.3 


0.14 


AzTEC/C143 d 


12.0 


J100149.80+015752.5 


0.066 


J100149. 77+015752.4 


1/6 


21.18+2.3 


0.43 










J100149.87+015731.8 


2/6 


12.45+3.9 


0.25 










J100149. 01+015735.4 


3/6 


12.25+3.0 


0.25 


AzTEC/C144 


16.1 


J100142. 55+020014.9 


0.012 


J100142. 54+020014. 8 


1/7 


56.60+6.6 


0.41 










J100142. 55+020019.8 


2/7 


45.93+21.7 


0.33 










J100142.13+020021.0 


3/7 


23.40+3.1 


0.17 


AzTEC/C145 


16.1 


J100031. 82+021243. 2 


0.023 


J100031. 84+021242. 9 


2/8 


58.68+2.1 


0.47 










J100031. 41+021232.5 


1/8 


60.08+8.0 


0.48 










J100031. 84+021242.9 


2/8 


58.68+2.1 


0.47 


AzTEC/C146 


16.1 


J095957.26+014818.0 


0.025 


J095957.22+014818.3 


4/7 


8.51+0.3 


0.11 










J095957.67+014803.7 


1/7 


31.37+2.9 


0.40 










J095957.31+014808.9 


2/7 


24.13+8.5 


0.30 


AzTEC/C147 d 


12.0 


J100107.80+015713.8 


0.015 


J100107.81+015712.5 


4/6 


3.93+0.8 


0.14 










J100107. 19+015723.6 


1/6 


10.40+3.7 


0.37 










J100107.40+015718.2 


3/6 


4.73+0.5 


0.17 


AzTEC/C153 d 


12.0 


J100212. 57+022431. 9 


0.013 


J100212. 52+022430. 8 


1/5 


43.47+8.2 


0.51 










J100212. 34+022433.1 


2/5 


37.11+8.1 


0.44 










J100212.42+022423.1 


3/5 


2.51+2.3 


0.03 


AzTEC/C155 


16.3 


J100054.59+020459.5 


0.095 


J100054.61+020459.2 


6/6 


0.35+0.1 


0.01 










J100054.48+020518.6 


1/6 


19.92+7.7 


0.48 










J100054.98+020512.1 


2/6 


14.86+1.8 


0.36 


AzTEC/C157 c 


16.4 


J100103.55+014810.7 


0.059 


J100103. 56+014810.8 


2/9 


27.02+0.8 


0.33 










J100104.76+014804.5 


1/9 


43.69+3.9 


0.53 










J100103. 56+014810.8 


2/9 


27.02+0.8 


0.33 


AzTEC/C159 


16.6 


J095930.40+015527.5 


0.075 


J095930. 40+015527.6 


3/8 


21.77+5.9 


0.17 










J095929. 58+015540.0 


1/8 


71.74+5.4 


0.57 










J095929. 80+015532.4 


2/8 


29.68+15.7 


0.23 


AzTEC/C160 


16.6 


J100023. 98+021749. 9 


0.005 


J100024. 02+021749. 7 


4/8 


3.73+0.2 


0.09 










J100023.82+021754.8 


1/8 


15.83+11.5 


0.39 










J100024.19+021756.8 


2/8 


13.70+4.2 


0.34 


AzTEC/C161 d 


12.0 


J095939. 11+023225. 7 


0.034 


J095939. 11+023224. 6 


1/10 


57.55+5.8 


0.55 










J095939.76+023213.8 


2/10 


18.67+9.6 


0.18 










J095940. 25+023218.2 


3/10 


11.19+8.7 


0.11 


AzTEC/C162 


16.6 


J100002. 21+021621. 9 


0.094 


J100002. 21+021621.9 


4/15 


3.63+0.5 


0.08 










J100001. 73+021604.0 


1/15 


19.81+10.6 


0.46 










J100002. 54+021609.8 


2/15 


6.16+0.4 


0.14 


AzTEC/C163 d 


12.2 


J100124.26+015847.9 


0.021 


J100124. 26+015847.8 


2/9 


7.16+1.6 


0.10 



Other IDs c 



AzTEC3 



COSLA-13 
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AzTEC ID a 


Rs 


Radio Coords b 


P 1 ,PH7 b 
i .4\jnz 


IRAC Coords b 


Rank b 


Likelihood 


Reliability 


Other IDs c 




(") 


(J2000) 


(J2000) 




Ratio b 


Factor b 












J100124.88+015849.0 


1/9 


52.68+7.2 


0.75 












J100124.65+015846.2 


3/9 


4.57+1.0 


0.06 




AzTEC/C164 d 


12.5 


J095935. 38+022236. 9 


0.009 


J095935. 31+022236. 2 


1/2 


52.28+9.0 


0.95 












J095935.35+022240.4 


2/2 


2.69+0.7 


0.05 




AzTEC/C167 d 


12.6 


J095922.83+020406.4 


0.126 


J095922. 81+020405. 8 


4/6 


0.80+0.4 


0.01 












J095923. 54+020351. 7 


1/6 


60.15+8.0 


0.84 












J095923.18+020356.1 


3/6 


1.72+0.4 


0.02 




AzTEC/C168 d 


12.6 


J100057.95+021311.4 


0.066 


J100057.98+021311.1 


4/7 


1.03+0.1 


0.01 












J100057.44+021308.3 


1/7 


65.05+10.5 


0.55 












J100057.97+021318.0 


3/7 


1.59+3.0 


0.01 




AzTEC/C171 


16.9 


J100047.53+020941.1 


0.002 


J100047.54+020941.1 


1/6 


23.02+3.4 


0.44 












J100048. 17+020930.5 


2/6 


10.27+1.6 


0.20 












J100047.04+020927.3 


3/6 


8.10+0.9 


0.15 




AzTEC/C173 d 


12.6 


J100118.61+015058.2 


0.053 


J100118. 64+015057.9 


6/8 


2.22+0.1 


0.03 












J100118.85+015048.1 


1/8 


59.96+5.9 


0.73 












J100117.97+015042.1 


3/8 


4.27+1.3 


0.05 




AzTEC/C174° 


17.0 


J100008. 27+015623. 9 


0.017 


J100008. 28+015624. 3 


1/8 


74.42+9.3 


0.58 












J100007.82+015614.4 


2/8 


32.73+5.1 


0.25 












J100007.64+015621.8 


3/8 


12.71+9.9 


0.10 




AzTEC/C176 


17.0 


J095957.11+021136.8 


0.013 


J095957.09+021136.4 


1/13 


49.45+2.7 


0.35 












J095957.26+021143.0 


2/13 


46.61+10.0 


0.33 












J095957.24+021127.0 


3/13 


16.85+14.0 


0.12 




AzTEC/C181 c 


17.2 


J095929. 55+021713.1 


0.047 












AzTEC/C186 


17.4 


J100105. 66+023238.1 


0.025 


J100105. 66+023238. 7 


1/8 


39.61+4.9 


0.73 












J100104.96+023242.0 


2/8 


7.65+3.1 


0.14 












J100104.99+023248.0 


3/8 


2.13+0.3 


0.04 





a AzTEC IDs correspond to the IDs from lDownes et al.| i2012ft for the GOODS fields and lAretxaga et alj i201l| ) for COSMOS. 
b Radio counterparts with P-statS$0.05 are bolded. IRAC counterparts that are ranked correctly (1st) and 
reliability factors greater than 0.8 are also bolded. 

c "AzTECxx" correspond to the IDs from [Younger et al.| j2007h and [Younger et al ] j2009l) for SMA observations, "COSLA-xx" correspond to IDs from 
ISmolcic et all <201Sft for PdBI observations of AzTEC sources in the COSMOS field, and "GOODS850-x" correspond to IDs from lBarger et al.| l |2012|) for 
SMA observations of SCUBA/AzTEC sources in GOODS-N. 

^Radio counterpart is identified in the extended 3er radio catalog (see Section |5.3t . 
e AzTEC source is cross- listed with Table [A~2l 
■^IRAC counterpart is not detected in all four bands. 

9 IRAC counterparts identified via the SMA (Younger et al. 2007, 2009). These include two radio-undetected SMGs (AzTEC/C4 and AzTEC/CIO) and one 
multiple system (AzTEC/C42), for which the SMA identified the correct counterpart and which is cross- listed with Tabic [A2l Two SMA sources, AzTEC7 
I Younger et al. 2007) and AzTECIO (Younger et al. 2009) do not match any SMGs in our AzTEC catalogs and so their IRAC fluxes were taken from their 
respective papers. 

'''The true radio counterpart for AzTEC/GN14 is undetected in our radio catalogs. The radio source listed here is unassociated with the submm emission 
iBarger et al.ll2012f ). 



Table A2: Multiple Radio Set: List of SMGs in GOODS-N, GOODS-S, and COSMOS that have multiple radio detections. 
The full table (108 rows) is available online. The IRAC counterpart to each radio detection is listed on the same line with 
its rank, likelihood ratio, and reliability factor. The information for an additional potential IRAC counterpart is often also 
listed. For COSMOS, radio detections which were found in the extended catalog (see Section 15.31 1 are indicated. Some 
AzTEC/COSMOS SMGs were single systems in the unextended radio catalog, but are multiple systems in the extended 
catalog. These sources were used in the original training set and are cross- listed with Table IAT1 



AzTEC ID a 


Rs 


Radio Coords b 


P-stat 14GHz b 


IRAC Coords b 


Rank b 


Likelihood 


Reliability 


Other IDs c 




(") 


(J2000) 


(J2000) 




Ratio b 


Factor 13 










GOODS-N 










AzTEC/GN7 


6.4 


J123711. 32+621330. 9 


0.028 


J123711. 34+621331.4 


3/4 


21.40+0.4 


0.17 








J123711. 99+621325. 5 


0.030 


J123711. 99+621326.0 


1/4 


48.23+2.4 


0.39 












J123711. 86+621333.8 


2/4 


45.87+3.4 


0.37 




AzTEC/GNll 


7.4 


J123635. 87+620707.8 


0.008 


J123635.85+620708.0 


1/8 


10.09+0.1 


0.37 








J123635. 87+620703. 9 


0.012 
















J123636. 32+620706. 8 


0.034 


J123636.37+620707.3 


3/8 


5.31+0.1 


0.20 




AzTEC/GN16 


7.9 


J123615. 94+621514. 3 


0.024 


J123615.82+621515.8 


1/5 


96.55+3.1 


0.51 


GOODS850-7 
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AzTEC ID a R s Radio Coords b P 1 4GHz b IRAC Coords b Rank b Likelihood Reliability 

(") (J2000) ' (J2000) Ratio 13 Factor b 







J123616. 10+621513. 8 


0.027 


J123616. 10+621514. d 


2/5 


74.93+1.2 


0.40 










J123616. 84+621514. 7 


3/5 


7.07+0.6 


0.04 


AzTEC/GN17 


8.0 


J123539.91+621430.8 


0.071 


J123539.96+621431.1 e 












J123539. 93+621441. 9 


0.076 


J123539. 94+621441. 17° 












J123541. 204-621438. 5 


0.034 


J123541. 25+621439. e 








AzTEC/GN18 


8.0 


J123741. 16+621220. 5 


0.020 


J123741. 16+621221. 3 


1/6 


95.07+3.6 


0.73 






J123741. 63+621223. 6 


0.047 


J123741. 66+621224.0 


2/6 


21.91+1.2 


0.17 










J123740.12+621220.0 


3/6 


4.80+12.6 


0.04 


AzTEC/GN24 


8.7 


J123608. 57+621435. 3 


0.058 


J123608. 60+621435. 7 


1/6 


33.28+1.5 


0.50 






J123608.64+621435.5 


0.056 


J123608.60+621435.7 d 


1/6 


33.28+1.5 


0.50 










J123607.66+621441.7 


2/6 


7.28+1.4 


0.11 


AzTEC/GN30 


9.1 


J123641. 85+621719.7 


0.068 














J123642. 18+621722.1 


0.049 


J123642. 19+621722. 7 


2/5 


7.05+0.1 


0.24 










J123642. 78+621712. 5 


1/5 


8.98+4.8 


0.30 


AzTEC/GN36 


9.5 


J123617. 39+621551. 2 


0.046 


J123617.44+621551.5 


4/6 


5.24+0.3 


0.12 






J123617. 54+621540. 7 


0.034 


J123617.55+621540.8 


1/6 


15.41+0.7 


0.35 






J123618. 32+621550. 5 


0.094 


J123618. 33+621550. 8 


3/6 


7.62+0.3 


0.17 










GOODS-S 








AzTEC/GS19 


10.3 


J033222. 71-274127.0 


0.080 


J033222. 71-274126.2 


3/7 


4.26+0.2 


0.13 






J033222. 92-274125. 3 


0.061 


J033222. 85-274125.1 


4/7 


4.03+0.1 


0.12 










J033223. 77-274131. 6 


1/7 


11.66+0.7 


0.34 


AzTEC/GS32 


13.7 


J033310. 14-275124.4 


0.195 


J033310. 12-275124.8 


3/4 


1.18+0.2 


0.20 






J033308.63-275134.5 


0.092 


J033308.62-275134.5 


1/4 


2.03+0.1 


0.34 










J033308. 78-275125. 2 


2/4 


1.50+0.4 


0.25 


AzTEC/GS35 


12.6 


J033227.19-274051.5 


0.004 


J033227.17-274051.7 d 


4/9 


27.92+1.1 


0.14 






J033226. 97-274049. 9 


0.009 


J033227.17-274051.7 


4/9 


27.92+1.1 


0.14 










J033226. 83-274056.1 


1/9 


50.66+8.3 


0.26 


AzTEC/GS37 


14.6 


J033255. 70-274623.8 


0.212 














J033257.09-274624.5 


0.150 


















J033256. 71-274606. 6 


1/5 


20.61+10.1 


0.81 










COSMOS 








AzTEC/C23 


11.0 


J100142. 36+021836.0 


0.013 


J100142. 34+021835. 2 


1/2 

1 


25.74+1.9 


0.90 






J100142. 76+021841. 5 


0.036 


J100142. 73+021841.1 


2/2 


2.62+0.2 


0.09 


AzTEC/C43 f 


12.5 


J100003. 12+020201. 5 


0.032 


J100003. 09+020201. 5 


5/8 


7.91+2.3 


0.05 






J100003.76+020217.6 8 


0.069 


J100003. 81+020216. 7 


7/8 


1.25+0.1 


0.01 










J100003. 40+020204. 4 


1/8 


65.08+7.8 


0.38 


AzTEC/C53 


12.8 


J100122. 47+021209. 7§ 


0.006 


J100122.46+021210.0 


4/8 


4.82+0.4 


0.05 






J100122. 07+021212. 5 s 


0.038 


J100122. 04+021213. 9 


2/8 

/ 


1.39+0.3 


0.02 










J100122. 61+021214. 9 


1/8 


48.66+14.8 


0.53 


AzTEC/C59 


13.2 


J100030. 13+023716. 7 


0.018 


J100030.14+023716.8 


1/5 


9.35+0.3 


0.51 






J100031. 11+023717.2 


0.066 


J100031. 10+023717.2 


4/5 


1.44+0.0 


0.08 










J100303. 32+014413. 7 


2/5 


3.29+1.0 


0.18 


AzTEC/C69 


13.5 


J100138. 00+020908. 7§ 


0.001 


J100138. 00+020910.0 


1/5 


31.40+4.3 


0.68 






J100137.29+020904.6 S 


0.076 


J100137.35+020904.3 


3/5 


4.53+1.2 


0.10 










J100138. 17+020913.3 


2/5 


8.27+1.6 


0.18 


AzTEC/C70 


13.5 


J100025. 47+020312. 6 


0.037 














J100025. 60+020316.1 


0.022 


J100025.60+020316.2 


3/4 


3.19+0.1 


0.09 










J100025. 94+020315. 2 


1/4 


27.53+12.9 


0.75 


AzTEC/C91 f 


14.5 


J100128. 20+022340. 8§ 


0.038 


J100128. 21+022340. 8 


3/8 


2.22+0.1 


0.07 






J100128. 48+022345. 2 


0.005 


J100128.48+022345.1 


1/8 


19.11+6.7 


0.57 










J100129.30+022348.2 


2/8 


7.33+0.5 


0.22 


AzTEC/ClOO 


14.7 


J095917.70+021028.8 


0.091 


J095917.67+021028.5 


4/7 


1.14+0.1 


0.02 






J095917.98+021043.7 


0.061 


J095917.96+021043.5 


2/7 


4.17+1.3 


0.06 










J095918. 70+021040. 9 


1/7 


61.71+10.5 


0.86 


AzTEC/CllO 


14.9 


J100108. 71+020036. 48 


0.001 


J100108. 72+020036. 6 


1/5 


35.27+2.6 


0.88 






J100108. 50+020029. I s 


0.025 


J100108.46+020030.8 


2/5 


2.15+0.6 


0.05 
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AzTEC ID a R s Radio Coords b P 1 4GHz b IRAC Coords b Rank b Likelihood Reliability 

(") (J2000) ■ - (J2000) Ratio 13 Factor b 











1100108 46+020030 8 


3/5 


O 1 K-Lf) ft 


0.05 


AzTEC/C124 


15.3 


J095946. 35+023602.1 


0.033 


J095946. 31+023602. 3 


3/6 


1.91+0.1 


0.09 






J095946. 68+023546. 8 


0.038 


T095946 63+023546 8 


1/6 


14 36+1 


0.68 










1095956 1 1 +01 4253 5 


2/6 


2 43+0 9 


0.12 


AzTEC/C140 f 


16.0 


J 100124. 96+015147. 3 


0.004 


J100125. 00+015147. 2 


2/5 


37.98+4.8 


0.35 






11 001 25 49+01 51 35 386 


0.071 


11 001 25 47+01 51 34 d 
















J100124. 60+015141. 3 


1/5 


61 43+9 8 


0.57 


AzTEC/C142 f 


16.0 


J100018. 82+020245.0 


0.070 


J100018. 83+020245.0 


6/7 


0.93+0.0 


0.03 






J100017. 46+020252. 6 g 


0.071 


J100017.47+020252.3 


2/7 


4.27+0.3 


0.14 










1100018 37+020245 1 


1/7 


18 97+14 7 


0.61 


AzTEC/C149 


16.2 


J100149. 82+022831. 9 


0.051 


J100149. 82+022831. 9 


4/9 


2.13+0.2 


0.06 






J 100150. 92+022826. 3 


0.026 


J100150. 92+022826. 4 


2/9 


8.46+4.0 


0.25 










1100149 96+022824 5 


1/9 


16 85+9 2 


0.49 


AzTEC/C150 


16.2 


J 100004. 81 +023045.1 


0.016 


J100004. 79+023045. 3 


1/5 


80.30+1.6 


0.91 






Tl 00005 43+023029 1 


0.089 


1100005 44+023029 2 


5/5 


28+0 


0.00 










1100005 74+023039 1 


2/5 


vJ . UO.L J. . VJ 


0.04 


AzTEC/C154 


16.3 


J 100035. 41 +023531. 8 


0.046 


1100035 41+023532 1 


4/10 


4.39+0.4 


0.09 






T100035 84+023543 8 


0.032 


1100035 85+023544 4 

•J -L KJ KIKJ ij Ul . O | 1 7 Zj O O a a . j: 


5 /10 


4 28+0 3 


0.08 










1100035 39+023541 


1/10 


19 16+1 6 


0.38 


AzTEC/C157 f 


16.4 


J100103. 55+014810. 7 


0.059 


J100103. 56+014810. 8 


2/9 


27.02+0.8 


0.33 






J 100104. 79+014805. g 


0.037 


1100104 76+014804 5 


1/9 


43 69+3 9 


0.53 










11 001 03 82+01 481 5 1 


3/9 


5.49+4.0 


0.07 


AzTEC/C158 


16.5 


J100013. 58+021230. 8 


0.104 


J100013. 58+021230. 8 


8/8 


0.06+0.0 


0.00 






T1 0001 5 28+021 240 fi 


0.085 


11 0001 5 28+021 240 2 

U iUTUi l^TUii 


3/8 


k 09+0 4 


0.11 










J100015. 01+021235. 2 


1/8 


29.52+24.0 


0.54 


AzTEC/C169 


16.9 


J100013. 36+023016. 2 


0.125 


J100013. 34+023016.0 


7/8 


0.28+0.0 


0.01 






T100014 21+023019 1 


0.051 


T100014 21+023019 


1 /8 

I/O 


27 68+0 6 


0.83 










J100013. 79+023019.1 


2/8 


1.82+1.4 


0.05 


AzTEC/C174 f 


17.0 


J 100008. 27+015623. 9 


0.017 


1100008 28+015624 3 


1/8 


74.42+9.3 


0.58 






T1 00007 50+01 5629 7S 


0.117 


1100007 54+015628 9 


8/8 


n 22+0 8 


0.00 










1100007 82+015614 4 


2/8 


32 73+5 1 


0.25 


AzTEC/C178 


17.1 


J 100141. 18+022002. 2 


0.010 


J100141. 15+022002.0 


1/10 


15 68+1 


0.46 






J100142. 30+022005. 3 


0.110 


J100142. 29+022005. 3 


5/10 


1.99+0.2 


0.06 










1100140 96+021957 7 


2/10 


q q«+i ^ q 


0.29 


AzTEC/C181 f 


17.2 


J095929. 55+021 713.1 


0.047 














T095929 47+021658 is 


0.119 


1095929 47+021658 5 








AzTEC/C183 


17.3 


J100226. 17+021227.1 


0.00 














J100226. 25+021229. 6 


0.003 


J100226.25+021228.8 


1/11 


7.92+0.5 


0.24 










J100225.70+021224.2 


2/11 


6.81+3.5 


0.20 


AzTEC/C184 


17.3 


J100015. 71+014445.6 


0.054 


J100015.70+014445.5 


5/10 


1.89+0.1 


0.04 






J100016.33+014437.4 


0.133 


J100016.34+014437.5 


7/10 


1.22+0.0 


0.03 






J100016. 77+014453. 9S 


0.044 


J100016.81+014453.18 


6/10 


1.72+0.1 


0.04 










J100015.42+014447.3 


1/10 


15.83+4.3 


0.37 


AzTEC/C187 


17.4 


J100152. 55+021954.3 


0.079 


J100152.57+021954.0 


2/13 


14.97+0.3 


0.28 






J100153. 34+021928.7 


0.115 


J100153.34+021928.7 


4/13 


2.66+0.0 


0.05 










J100153.46+021939.6 


1/13 


20.67+5.0 


0.39 



"AzTEC IDs correspond to the IDs from lDownes et al.| j2012ft for the GOODS fields and lAretxaga et all j201lj ) for COSMOS. 
b Radio detections with P-statsCO.05 are bolded. IRAC detections with reliability factors greater than 0.8 are also boldcd. 

c "AzTECxx" correspond to the IDs from Younger ct al. (2(307) and [Younger et al for SMA observations, "COSLA-xx" correspond to IDs from 

ISmolcic et al.| j2012f) for PdBI observations of AzTEC sources in the COSMOS field, and "GOODS850-x" correspond to IDs from lBarger et al.| j2012|) for 
SMA observations of SCUBA/AzTEC sources in GOODS-N. 
^IRAC detection is within l" of both radio detections. 
e IRAC counterpart is not detected in all four bands. 

AzTEC source is cross- listed with Table [All 
9 Radio counterpart is identified in the extended 3tr radio catalog (sec Scction |5.3| . 
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Table A3: Radio-Undetected SMG Set: list of SMGs without radio detections. The full table (321 rows) is available online. 
For each AzTEC source, the likelihood ratios and reliability factors for the top three ranked potential IRAC counterparts 
are listed. IRAC detections with R ^ 0.8 are bolded. 



AzTEC ID a 


Rs 


IRAC Coords b 


Likelihood 


Reliabili 




(") 


(J2000) 


Ratio b 


Factor 






GOODS-N 




AzTEC/GN2 


4.8 


J123632. 01+621711.1 

~ 


15.26+2.1 


0.99 

— 


AzTEC/GN8 


6.6 


J123645. 85+621447.3 


12.23+0.4 


_ 
0.70 






J123645. 77+621442. 6 


5.05+0.6 


0.29 


AzTEC/GNIO 


7.3 


J123627. 54+621218. 4 


59.47+3.7 


— 
0.63 






J123627. 01+621217. 7 


19.76+14.6 


0.21 






J123627. 18+621212. 9 


6.10+11.9 


0.06 


AzTEC/GN13 


7.7 


J123554. 28+621343. 8 


80.05+3.8 


0.64 






J123553.57+621341.3 


25.17+5.0 


0.20 






J 123554. 54+621347. 7 


12.62+2.5 


0.10 


AzTEC/GN21 


8.3 


J123800.84+621612.2 


23.37+0.8 


0.57 






J123800. 41+621607. 4 


7.34+8.0 


0.18 






J123801. 77+621611. 6 


3.88+1.8 


0.10 


AzTEC/GN22 


8.5 


J123649. 60+621212. 9 


8.64+8.8 


0.28 






J123650. 39+621209. 9 


8.53+5.5 


0.28 






J123649. 37+621214. 8 


6.95+9.2 


0.23 


AzTEC/GN27 


8.8 


J123720.07+621223.2 


70.18+2.6 


0.49 






J123719. 61+621221. 3 


33.89+1.0 


0.24 






J123718.95+621217.7 


29.71+1.9 


0.21 


AzTEC/GN28 


8.9 


J123644.03+621938.8 


71.77+6.2 


0.68 






J123642.99+621931.0 


23.21+12.2 


0.22 






J123643.98+621930.8 


6.74+6.3 


0.06 


AzTEC/GN29 


9.1 








Az 1 rLO / <jlNol 


9.1 


J123622.07+621616.2 


93.97+3.9 


0.59 






J123622.49+621609.4 


42.17+34.8 


0.27 






J123621.82+621609.0 


15.37+14.7 


0.10 


A 'TT~ f ' //^MOO 

AzlrjO/OlNdi 


9.1 


J123716. 87+621401.1 


79.87+16.4 


0.87 






J123717.16+621352.5 


5.44+0.1 


0.06 






J123717.28+621356.8 


3.19+0.2 


0.03 






GOODS-S 




AzTEC/GS2 


5.3 


J033218.28-275224.5 


11.01+0.6 


0.69 






J033218. 84-275223.3 


4.71+0.6 


0.30 


AzTEC/GS3 


6.1 


J033247.96-275416.3 


27.53+1.9 


— 
0.55 






J033247.69-275423.4 


22.07+4.6 


0.44 


AzTEC/GS4 


6.1 


J033248. 93-274251. 3 


92.64+6.6 


1.00 

- 


AzTEC/GS6 


8.0 


J033225.77-275220.0 


1.61+0.1 


0.59 






J033226.25-275222.0 


0.52+0.0 


0.19 






J033225.52-275226.7 


0.39+2.7 


0.14 


AzTEC/GSll 


9.2 


J033215. 29-275038. 8 


2.37+0.1 


0.81 






J033215.37-275044.7 


0.32+0.0 


0.11 


AzTEC/GS12 


9.3 


J033229.29-275619.0 


64.80+4.9 


0.78 






J033228.91-275614.1 


18.39+27.5 


0.22 


AzTEC/GS13 


9.0 


J033211. 92-274615.2 


30.42+3.6 


0.72 






J033211. 56-274613.0 


4.74+0.3 


0.11 



Other IDs 



GOODS850-11 



GOODS850-2 
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AzTEC ID a R s IRAC Coords b Likelihood Reliability Other IDs c 

(") (J2000) Ratio b Factor b 









91 974fi90 & 


9 qi _|_n 1 


07 


A7TFC/OSI4 


9.2 


1033234 


75-27521 7 4 


q fiS+6 5 

y . uo_i_u. o 


0.39 






T033935 


01 97^91 Q ^ 


Qfi-I-I 1 
U.OUX1. 1 


9fi 
U.zu 






T033934 


9k 97^99^ Q 


4.31i0.4 


0.18 


AzTF.C/nSI 5 

ili -L IJVj / VJkJ J- U 


9.4 


T0331 51 


S7-274601 1 


1 9 41 -h5 

_L Z . rr _1 1_ O • U 


0.59 






T0331 50 

J UOO-LOU 


70 974^0^ 4 


5.71il.7 


0.27 






T0331 50 


Qfi 974^00 7 


O fiQ-UO 4 


13 

u.io 




10.0 




QQ-2744fln 8 

ij \ J Z« f 'l^UUiO 


47 5«-U9 7 

^± J iUO 1 Zrf • I 


0.91 






T033937 


974407 


1.82±0.1 


0.03 






T033937 


974404 9 


1.57i0.1 


0.03 




10.0 


T033999 


1 97481 1 


1 Q 44-I-0 R 
± y . . u 


0.55 






T033999 


^9 97481 ^ ^ 


i q 77-1-0 3 

LO.il ^U.O 


3Q 

u.oy 






T033999 


1 ^ 974S99 


1 4Q-I-3 9 


04 
U.U1 


A7TFC/GS18 

AZi _L lJU / V_XlO lU 


9.6 


T033243 


^9-9746^ 

. OZi Z f . \J 


qq 73+4 4 


0.67 






T033944 


00 974fi' : t^ 


QQ 1 QJ-1 9 

oy . ±o^-L .z 


98 
u.zo 






T03394 3 


44 974^*^4 ^ 


c; 79-1-0 1 

O. 1 Z^U. 1 


04 

U.U1 


A7TFC/OS22 


12.0 


T03321 2 


1 5-974949 Q 


1 9 QS-hO (S 

_L Z . .VI ' 1 1 1 , U 


0.39 






T03391 9 


54 974^08 


1 1 94-1-0 7 


34 






T03391 9 


5fi 974959 


c; 11 _|_r; r: 


15 




11.9 


T033991 


1 ^ 975fi9fi 5 


75 31 -1-3 9 
■ o . o±^o . z 


0.75 






T033991 


(SO 975fi9 ,: i ^ 
■ UU z / ouzo.o 


1 1 ni _|_n o 

± _L . U-L ZCU . Z 


0.11 






T033991 


0Q 975fi90 7 


4 Q9-U1 c; 


05 

u.uo 


iAZj _L 1_jV^ / VjxQZI 


11.8 


TO" : i' : 49 i: izL 

J UOOZO^i 


9^ 974Q^Q S 
. zo z / lyos.o 


7.27±0.4 


0.65 






T033935 


1 1 974Q4S 4 


9 QQ-l-1 


9fi 
u.zu 






T033934 


fiQ 97/1QQQ n 


n ^Q-l-O 1 


05 
u.uo 


A7TFP /rm^fi 


11.6 


T03391 


■^0 974343 3 


J. y . o / . u 


0.49 






T03391 


95 974330 5 
zo-z / IOOU.U 


4 QQ-l-1 5 
i.yyx iu.u 


12 






T03391 5 


fi5 974330 ^ 
UO-z / IOOVJ.O 


4 45-1-5 4 


11 


A?TFP /CN97 

-TlZr _L JZj / VjOZ / 


12.2 


T033941 


fi4 9741 59 1 

. Urt Z 1 '-t _L O Z . J. 


n Q74-0 


0.19 






T033949 


33 974904 1 


n eq-UO 1 

U.OOxU. 1 


1 s 
u. io 






T033949 


7Q 9741 44 & 
. i y-z / L t± L t i -t.o 


n 70-1-0 3 


15 


/\Zr _L JZj O / VjlJiO 


12 1 


Tfl^**949 

•J UOOi'l/ii 






0.89 






T033949 


Q0 97591 Q Q 
. yu z / oz ± y . y 


3 35-1-9 8 
o . oox^. o 


0.04 






T033949 


74 975905 
. I i-z / ozuo.u 


9 7Q4-0 5 
Z. ( o^u.o 


03 

u.uo 




12.1 


T0331 5Q 

U UOO IUj 


31-974500 5 


1 5 72+fi ft 


0.58 






J033159 


.71-274459.8 


7.46±3.0 


0.28 






J033158. 


.68-274500.3 


3.46±0.3 


0.13 




13.5 


T033990 


fifi 974934 5 


9 29+0.2 


0.29 






T033991 


9fi 974945 Q 
,zu _ z / izu.y 


f. c;4+c; 1 

0.01X0. 1 


90 
u.zu 






T033991 


51 974941 & 


5 &7-I-0 1 
O.O / IU. 1 


18 
u. io 


A7TFP /RSffi 


13.6 


•J Ui)Odl4i 


JO Za l O1JJ.1J.O 


q-i ccxq i 


0.91 






T03391 3 


9Q 97551 Q Q 
zy-z / ooiy.y 




04 
U.U1 






T03391 4 


1 1 975594 Q 


9 S4+S 7 


03 
U.UO 


A7TFP /d^dO 


14.6 


T033901 

J UOOZU-L 


1 5 974fi3fi 


3.74zb0.6 


0.64 






T033900 

J uoozuu 


(S3 974fi33 1 

.UO Z / 1UOO. ± 


i Q«+n 3 


0.32 


A7TFC/fiS4S 

AZj _L Ljv / VJXiJrtO 


12.9 


T03321 fi 

U VjvJOii J-U 


53-97594fi 5 


X 9Q+0 3 


0.32 






T03391 


04 975954 5 


7 77+5 


30 

u.ou 






T03391 5 


74 975943 4 


e; QQ+n 9 

O . OOIU . z 


0.23 


A7TFC/GS49 


13.7 


T033299 


03-274003 6 


p; 3S+0 7 


0.55 






T033991 


9Q 97400Q K 
.zy-z / luuy.o 


1.70i0.1 


0.17 






T033991 


1 97401 Q 4 
. iu _ z / ivjiy .i 


1 34+0 1 


14 


A7TFC/OS5O 

AZj _L Lju / V_XkJtJU 


13.7 




01-974Q39 


3.87±2.7 


0.34 






JUoozoo. 


,4y-z / 4yz0. 1 


o.zuitl. r 


U.zy 






J033235 


.36-274920.7 


2.73±0.6 


0.24 


AzTEC/GS52 


15.4 


J033243 


.59-275011.6 


14.71±10.4 


0.42 






J033244 


.45-275014.9 


10.92±17.7 


0.31 






J033243 


.90-275017.3 


2.00±8.3 


0.06 


AzTEC/GS53 


15.7 


J033204 


.01-274453.6 


44.02±11.5 


0.64 






J033204 


.87-274451.7 


11.72±2.7 


0.17 






J033205 


.12-274501.4 


3.79±0.9 


0.05 


AzTEC/GS55 


16.0 


J033225 


.74-275818.2 


1.04±0.1 


0.37 






J033224 


.93-275825.7 


0.76±0.4 


0.27 
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AzTEC ID a 


R s 


IRAC Coords b 


Likelihood 


Reliabil 




(") 


(J2000) 


Ratio b 


Factoi 






J033225. 11-275818.2 


0.35+0.0 


0.13 






COSMOS 




AzTEC/C5 


6.9 


J095942. 85+022938. 2 

- 


87.20+14.5 

— 


1.00 


AzTEC/C7 


7.7 


- 

J100015.89+022943.9 


- 

8.87+12.4 


0.49 






J100015. 74+022939.6 


7.44+0.7 


0.41 






J100015. 48+021550.7 


0.93+0.2 


0.05 


AzTEC/Cll 


8.4 


J100141. 19+020357.1 
— 


0.03+19.2 
— 


0.11 


AzTEC/C19 


10.9 


- 

J095949. 91+015329. 9 

- 


— 

11.76+7.6 

- 


0.98 


AzTEC/C21 


11.0 


- 

J095921. 33+022238.5 


- 

3.99+0.2 


0.31 






J095921. 45+022232.1 


2.70+1.5 


0.21 






J095921. 06+022234.9 


2.00+1.0 


0.16 


AzTEC/C26 


11.2 


J100132.09+023215.9 


10.16+1.3 


0.52 






J100132.38+023215.0 


7.88+7.9 


0.41 






J100131. 85+023206.7 


0.58+0.0 


0.03 


AzTEC/C27 


11.2 


J095936. 48+020656.7 


5.18+2.9 


0.52 






J095936. 48+020656.7 


2.32+1.7 


0.23 






J095936. 78+020655.4 


2.09+3.4 


0.21 


AzTEC/C30 


11.4 


J100025. 42+022610.5 


11.82+2.4 


0.50 






J100025. 25+022606.0 


5.47+0.8 


0.23 






J100024.75+022604.7 


3.53+0.7 


0.15 


AzTEC/C34 


11.9 


J100007.54+021147.0 


59.07+9.9 


0.43 






J100007.45+021152.0 


38.02+4.2 


0.28 






J100008. 18+021153.6 


26.40+9.3 


0.19 


AzTEC/C39 


12.1 


J100126. 52+020006.1 


52.68+5.8 


0.81 






J100126.83+020001.7 


4.10+1.7 


0.06 






J100127.19+020000.8 


3.79+0.3 


0.06 


AzTEC/C49 


12.6 


J100131.80+015412.1 


12.17+3.4 


0.36 






J100131.33+015406.1 


7.80+1.3 


0.23 






J100132.03+015357.2 


6.01+1.7 


0.18 


AzTEC/C51 


12.7 


J100040.59+015923.5 


39.61+13.3 


0.57 






J100039.99+015924.6 


6.21+1.1 


0.09 






J100040.40+015927.5 


5.30+6.1 


0.08 


AzTEC/C54 


12.9 


J100126.01+015751.4 


62.24+10.5 


0.51 






J100126.22+015742.9 


38.44+25.5 


0.31 






J100126.51+015743.7 


17.08+2.8 


0.14 


AzTEC/C56 


13.0 


J095905.03+022156.8 


48.69+4.2 


0.69 






J095904.92+022149.7 


13.62+13.2 


0.19 






J095905. 26+022153. 3 


4.64+0.4 


0.07 


AzTEC/C60 


13.2 


J100128.39+022127.2 


55.19+18.3 


0.47 






J100128.37+022134.3 


27.20+21.1 


0.23 






J100128.19+022121.6 


21.92+14.2 


0.19 


AzTEC/C62 


13.3 


J100100.38+023802.8 


78.86+8.9 


0.72 






J100100.46+023805.9 


23.90+4.3 


0.22 






J100059.89+023756.8 


6.81+2.4 


0.06 


AzTEC/C63 


13.3 


J095920.05+023110.2 


6.58+2.5 


0.79 






J095921.07+023123.3 


0.64+0.3 


0.08 






J095920.26+023116.1 


0.54+3.5 


0.06 


AzTEC/C64 


13.4 


J100139. 32+022341.4 


73.90+14.5 


0.83 






J100139.54+022339.2 


4.27+0.8 


0.05 






J100139.88+022347.7 


2.43+14.2 


0.03 


AzTEC/C67 


13.5 


J100118.58+020933.0 


16.85+2.8 


0.76 






J100118.40+020945.0 


2.56+0.8 


0.11 






J100118.79+020950.3 


2.17+1.8 


0.10 


AzTEC/C68 


13.5 


J100121. 65+020728. 3 


4.90+1.0 


0.84 



Other IDs c 



AzTECl 
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AzTEC ID a R s IRAC Coords b Likelihood Reliability Other IDs c 

(") (J2000) Ratio b Factor b 



AzTEC/C73 13.7 

AzTEC/C75 13.8 

AzTEC/C76 13.8 

AzTEC/C78 13.9 

AzTEC/C80 13.9 

AzTEC/C81 13.9 

AzTEC/C83 14.2 

AzTEC/C89 14.4 

AzTEC/C90 14.4 

AzTEC/C92 14.5 

AzTEC/C94 14.6 

AzTEC/C96 14.7 

AzTEC/C99 14.7 

AzTEC/ClOl 14.8 

AzTEC/C102 14.8 

AzTEC/C104 14.8 

AzTEC/C106 14.9 

AzTEC/C107 14.9 

AzTEC/C108 14.9 

AzTEC/C114 15.1 

AzTEC/C117 15.2 



J100123.06+020730.5 
J100121.63+020731.5 
J100103.34+022845.7 
J100103.75+022853.3 
J100103.89+022907.6 





13+020059.7 


1100052 


.61+020112. 


3 


11 00053 

<J -1 UUUiJ'J ■ 


.24+020105. 


1 


11 0001 2 


.92+021211. 


6 


11 0001 3 


.35+021200. 


5 


11 0001 3 


.79+021203. 


1 


T0Q^Q03 


.21+015938. 


•1 


T0Q5Q02 


.24+015949. 


6 


T0Q<^Q09 
j uyoyu.^ . 


.36+015943. 


6 


ti ooos^ 


.40+022555. 


3 


ti 00033 

J 1UUUOO. 


.35+022601. 


7 


Ti 00032 

<J -i- \JVI\J' > ■ 


.55+022546. 


6 


TI 0000^ 

J -LUUUUO ■ 


.47+015241. 





TI OOOOfi 

J 1UUUUU. 


.11+015239. 


7 


TI OOOOfi 

J 1UUUUU. 


.06+015233. 


2 


TI 0022Q 


93+021416.9 


TI 00930 


.40+021409. 


6 


TI 0022Q 


.73+021419. 


8 


TI 001 97 


.34+021339. 


3 


Ti 001 27 


.69+021339. 


1 


TI 001 97 


.18+021332. 


9 


TI 001 S 1 ! 


.21+021649. 


•1 


TI 001 3^ 


.79+021658. 


8 


TI 001 Sfi 


.08+021652. 


8 


TI 001 40 


.48+023010. 


5 


TI 001 40 


.06+023019. 


5 


TI 001 3Q 

J 1UU 103. 


.16+023018. 


9 


T0Q^Q^7 
juyoyo i . 


.40+021715. 


3 


J095956. 


.80+021728. 


6 


J095957. 


.54+021729. 


8 


J100108. 


.93+015150. 


6 


J100109. 


.21+015155. 


9 


J100108. 


.23+015145. 


2 


J100006. 


.99+015957. 


6 


J100006. 


.81+015959. 


9 


J100007. 


.17+020008. 


1 


J095945. 


.34+023016. 


6 


J095945. 


.35+023027. 


1 


J095945. 


.21+023021. 


•1 


J095948. 


.12+023160. 





J095948. 


.91+023205. 


7 


J095949. 


.22+023145. 


7 



J095940.40+015336.0 

J095941.30+015336.7 
J095941. 30+015336.7 
J100007.05+023848.6 
J100006.16+023831.1 
J100005.96+023831.3 
J095939.65+022242.0 
J095939.87+022233.5 
J095939.49+022247.2 
J100115.67+023611.2 
J100116.29+023607.4 
J100116.35+023611.3 
J100024. 93+022004. 8 
J100023.62+022002.5 

J095926.59+022018.7 



A(\A-A A 


OS 

u.uo 


0.24zbl.3 


0.04 


^ 00-1-9 


44 


A 39-1-0 9 


0.38 


1.58±0.7 


0.14 




0.98 


1 93-1-0 1 


01 


f) OQ-LO 
U . ZOIU . u 


0.00 


19 33+3 9 


0.67 


a qi _|_n q 

U . O-L XlU . O 


0.13 


o 37-1-0 ^ 


0^ 

u.uo 


3 7^-1-0 


3Q 

u.oy 


2.51±2.4 


0.26 


9 4^-l-S 1 


9fi 


U-L .U-LXIU.O 


^7 


OO. 1 o^o.o 


3fi 
u.ou 


4.72±2.5 


0.04 


ZU.O-LXI^.U 


fil 


1 3 fi4-M 7 


31 


9 70-l-Q S 
z. ; ux:y .o 


Ofi 

u.uu 




0.83 


/I ni -I-90 Q 

i.ui^zu.y 


07 
u.u / 


3.15±1.3 


0.05 


qq 90-1-1 Q 


(\7 

u.u / 


e; 99-1-0 9 


OQ 

u.uy 


/I QQ-l-1 


OS 
u.uo 


A(\ QQ.-U/1 Q 




u.oy 


93 c.fi-i-o n 

zo.uuxo.u 


34 


Q (SO-r-0 S. 


0.05 


Qp. QQ-l-9 ft 


^ 
u.oo 


1 ^ 31 -l-Q 


93 


97-1-9 


OQ 
u.uy 


16 73±1 7 


0.74 


9 9Q_|_7 7 


10 
U. 1U 


i 77-1-0 9 


OS 
u.uo 


Q 37-U( : l 9 
y .o / zcu.z 


^4 

U.O'i 


4 81il 4 


9S 
u.zo 


1 QQ-UO 9 


11 
U. 1 -L 


it; n n 
lo.uix^t.y 


3Q 

u.oy 


i i cc;_|_7 c 
1 J..UJX / .U 


9Q 
u.zy 


Q Q3_|_7 Q 

OiOOi / .y 


22 


i c; 7^-1-7 
li). 1 ox: / .u 


4^ 


it; i c^-l-fi q 
id. idxu.y 


44 


3 (\AA-C\ 9 


10 
U. 1U 


K q«_|_-i 3 

. 3UI1 . 


0.34 


e; c:i _|_n ^ 


39 
u.oz 


1.94zbl.8 


0.11 


1 R 04.-1-1 fi 1 


0.81 


1.71i6.5 


0.09 


&7-I-0 9 
u.o / x;u.z 


04 

u.u^ 


7.40±3.2 


0.30 


R y1fi_l_9 tr 


n 97 


4.19±2.7 


0.17 


19.04±6.4 


0.50 


13.96±6.7 


0.37 


2.25±0.1 


0.06 


63.46±7.1 


0.62 


20.28±0.3 


0.20 


6.23±0.5 


0.06 


1.27±0.4 


0.82 


0.00±1.2 


0.00 



19.02+1.7 0.48 



© 2013 RAS, MNRAS OOO.ITHSol 



Submm/mm Galaxy Counterpart Identification 29 



AzTEC ID a R s IRAC Coords b Likelihood Reliability Other IDs c 

(") (J2000) Ratio b Factor b 
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30 Alberts et al. 



AzTEC ID a R s IRAC Coords b Likelihood Reliability Other IDs c 

(") (J2000) Ratio b Factor b 
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a AzTEC IDs correspond to the IDs from lDownes et alj |2012|) for the GOODS fields and lAretxaga et alj j201lj ) for COSMOS. 
b IRAC detections with reliability factors greater than 0.8 are boldcd. 

c "AzTECxx" correspond to the IDs from Younger ct al. ( 2(3(37) and Young er et al for SMA observations, "COSLA-xx" correspond to IDs from 

ISmolcic et alj <201Sfl for PdBI observations of AzTEC sources in the COSMOS field, and "GOODS850-X" correspond to IDs from lBarger et alj l2012f ) for 
SMA observations of SCUBA/AzTEC sources in GOODS-N. 
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